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Abstract 

In this paper we study the stability of the unique continuation in the case 
of the wave equation with variable coefficients independent of time. We prove 
a logarithmic estimate in a arbitrary domain of M n , where all the parameters 
are calculated explicitly in terms of the C 1 -norm of the coefficients and on the 
other geometric properties of the problem. We use the Carleman-type estimate 
proved by Tataru in 1995 and an iteration of the local stability. We apply the 
result to the case of a wave equation with data on a cylinder an we get a stable 
estimate for any positive time, also after the first conjugate point associated 
with the geodesics of the metric of the variable coefficients. 
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1 Introduction 

We consider the wave operator in M n+1 , 

n n 

P(y, D ) = -£>o + E 9 jk (x)DjD k + E h j {x)Dj + q(x), (1.1) 

j,k= 1 j=l 

where i/ = (l,i)6lx M n are the time-space variables, Dq = —idt, Dj = —id Xj . The 
coefficients g jk G C 1 (M n ) are real and independent of time, and [g ]k ] is a symmetric 
positive-definite matrix. The coefficients h\q G C 0 (M n ) are complex valued and 
independent of time. 

An operator P(y, D) is said to have the unique continuation property if for any 
solution u to Pu = 0 in a connected open set 12 C R” +1 and vanishing on an open 
subset B C 12, it follows that u vanishes in 12. 

In the paper [28] Tataru proved for the first time the unique continuation property for 
fll.ip across every non-characteristic C 2 -hypersurface with no limitation to the normal 
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direction. The key point of these results is a Carleman-type estimate involving an 
exponential pseudo-differential operator. 

Much is known about the consequences of the general unique continuation property 
for the corresponding Cauchy problem. Actually the unique continuation property 
has proved to be instructive in many areas of mathematics, e.g. in studying the 
uniqueness for linear and nonlinear PDEs together with their blow up or traveling 
wave solutions HQ, in studying the Anderson localization [7], in control theory to get 
controllability results [30] 0Tj, in inverse problems to obtain uniqueness and stability 
estimates [18] . In particular Tataru’s result [28] is crucial for the development of the 
Boundary Control method (see [5] for pioneering work and H3 for detailed exposition 
of the further developments). 

Concerning the continuous dependence of the unique continuation property, that is 
its stability, less results are available. The elliptic and the parabolic cases have been 
studied in several settings by using either Carleman estimates or some versions of the 
three ball theorem (see [T], for a review of the results). 

To our knowledge the hyperbolic case like (II.ip is still open for arbitrary domains 
and arbitrary matrix valued coefficients g lk (x), while there exist results for particular 
coefficients or domains (see [21] [32]). This is maybe related to the difficulty of using 
the standard Carleman estimates for hyperbolic operators in order to prove the unique 
continuation close to the characteristic directions, that is the reason why Tataru’s 
work was so important in this field. 

The aim of the present work is then to prove a global stability estimate for the unique 
continuation of the operator P(y, D). 

In a previous work [5] we proved this property for the local case. Namely, given S = 
{y G = 0} a C' 2,p -smooth oriented hypersurface, which is non-characteristic 

in hi, for some fixed p G (0,1), we assume that u G H l {yt) is supported in {y,tp(y) < 
0} D fi, and P(y,D)u G L 2 (Q). Then, for each y 0 G S, with i/)'(yo) 7^ 0, we find R,r 
with R > 2r > 0 such that the following stability estimate holds: 


u\\L 2 {B(y 0 ,r)) < hill 


u\\m(B(y 0 ,2R)) 


111 (l + 


H“llg 1 (Bfa n ,2fl)) 

W Pu Wl 2 (B( Vo ,2R)) 


Here B(y 0 ,r) is a ball in M n+1 of radius r > 0 centered in y 0 and B(y 0 ,r) C 
B(y 0 ,2R ) C fi. The radii r and R and the coefficient Cm have been explicitly 
calculated with dependency on the geometric parameters and on the function in 

i. 

In this work we use the previous local stability inequality to prove a similar logarith¬ 
mic estimate for quite general domains of K n+1 . 

Moreover we propose a procedure to calculate all the constants involved, dependent 
on the norms of i / j , the coefficients in fll.ljl . the properties of the domains and the 
smooth localizers. The procedure is described in Proposition 12.41 and Appendix A. 
Concerning the proof, in the unpublished manuscript [27], Tataru suggested the pos¬ 
sibility of obtaining a log-stability result, by splitting the estimate for high and low 
temporal frequencies and by using Gevrey-class localizers to improve the estimates 
of u for low temporal frequencies. Here and in [B] we have advanced that idea, by 
employing tools of subharmonic functions and proper choice of the localizers in the 
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iterating procedure, together with the explicit computations of the uniform radii r, R 
and the time frequencies used in the iteration. Of fundamental importance is the 
calculation of the positive lower bound of the radius r, without which the iterative 
procedure could stop before covering the desired subdomain of fh 
The technique used consists in iterating the local stability result, but considering 
the low temporal frequencies separately from the high temporal frequencies. The 
advantage is that one can avoid the usual (In In ...| In ||Pw|| lAY 9 iterated estimate 
(for ||tt||#i = 1 and ||P «|| L 2 << 1, 6 e (0,1)) and get a (| In ||Pu|| i 2 |)" 61 results. As a 
consequence one obtains a stable control of the solution u inside Q. for any positive 
time. Moreover, we can come as close as we wish to the optimal time of the control 
T opt , i.e. the time to reach the uniqueness in (see Corollary 13.41 as example of 
computation). The importance of this issue has also been underlined in [24], who 
worked with FBI transform technique to get a log-stability estimate for large times. 
Hoermander in [13] proved an upper bound of the type y/27/23T opt . The issue of 
reaching T opt for (11.111 has been solved in (28], see also [14] [25]. Here we can derive 
the stable determination of it. 

Like in the elliptic case, many possible applications can be derived out of it. In 
particular we plan to use the inequalities in Theorem I1.HI1.2I to obtain an explicit 
modulus of continuity for the inverse problem for the wave operator on manifolds. 
This would improve the existing inverse stability results for Riemannian manifolds, 
which are currently based either on compactness-type arguments, see in , or on 
very strong geometrical conditions for the coefficients, e.g. in [101 [EH, 2f|. Here is 
important to be able to relate the explicit estimates with some geometric invariant 
of the manifold (Ricci curvature, sectional curvature, diameter, etc.). 

As application, in section |3] we apply Theorem 1 1.21 to the case of an arbitrary domain 
of influence in M n+1 . This is a special case of manifold, once one considers g jk as the 
inverse of the metric tensor. We start with a time-cylinder where the wave solution 
vanishes (or has small data) and we get the stability in any compact subset of the 
associated domain of influence at time T. The control of solution in a stable way in 
the domain of influence can have numerous important applications in inverse prob¬ 
lems and in in control theory. Here we consider also the case in which the ray held 
has also singularities, i.e. behind the corresponding cut-locus. This means that in 
principle we are able to deal with manifolds that possesses conjugate points, trapped 
rays and other singularities of geodesics. Thus, we remove the usual non-trapping 
conditions used in the Carleman estimates. 

The paper is organized as follows: in Section [2] we prove Theorems 11.11 and 11.21 in 
Section |3] we present the application to the case of a domain of influence of the wave 
solution vanishing in a small cylinder. In Appendix A we present the table with the 
estimates for the parameters used in Sec. [2] and we study the uniform estimates for 
the distance function d g and the related function 0 defined in Sec. |3] 

We first introduce some assumptions. 

Assumption A1 Let be a connected open subset of R x R n . Let P(y, D) be the 
wave operator (II.ip . with g' jk {x) e C 1 (fl), h j , q e C'°(fl). We assume that u <G H l {VL) 
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and that P(y,D)u E L 2 (fl). Assume that there is a function if> E C 2,p (fl), for some 
p E ( 0 , 1)5 such that in a domain fl 0 C fl one has p{y,fi>'(y)) 7 ^ 0 and ^'(v) 7 ^ 0) 
where p(y,f) = —+ 9^ k { x )ijik is the principal symbol of P. 

Assume that there exist values fimin < max and a connected nonempty set T C flo 
such that: supp (u) fl T = 0; and 0 7 - {y E flo > 'fi max } C T (which implies 
that fl 0 contains a subset T where u vanishes, and that the value fi> m ax is obtained 
for points inside the domain fl 0 ). 

Assume that fi mm is such that the open set fl a = {y E fl 0 \T : if min < fi(y) < if ma x } 
is nonempty, connected and satisfies dist(<9flo, fl a ) > 0. 

See remark [2781 for comments about the construction. 

Assumption A2 We define A(Dq) to be a pseudo-differential operator with symbol 
a(£o) ; 0 < a < 1, where a E Cq°(R) is a smooth localizer supported in |£ 0 | < 2 , equal 
to one in |£ 0 | < 1- Furthermore let a E (^“(K) for a hxed a E (0,1). Here G ] f° is 
the set of Gevrey functions of class 1/a with compact support, defined in [151 23]- 
We also define the smooth localizer b(y), supported in \y\ < 2, 0 < b < 1 and equal 
to one in \y\ < 1. 

The main results of the paper are the following Theorems 11.11 and 11.21 together 
with their application in Section [3] Theorem 13.31 


Theorem 1.1. Under the conditions of Assumption A1-A2, define the open set fli = 
fl 0 \Y containing fl a . Then for every 0 < 9 < 1 we have 


rt|U 2 (n a ) < Ci 60 




(In (l + 


IMIgqni) \\Q' 


Moreover, for any m E (0,1] we get 
||w||ffi-m(n 0 ) < cff 0 


Ml/wni) 


In (l + 


“ ffl 




ll- Pu lh 2 (n 1 ) 


)) 


m6 


The constant Ci 60 is calculated in the proof. 


The dependency of the constant C 160 from the geometric parameters of the prob¬ 
lems and from ip and 6 is described in Proposition 12.41 


Assumption A3 Let fi be a connected open subset of K x M n . Let P(y, D) be 
the wave operator (11.11) . with g jk (x) E C 1 (0), h ?, q E C°(fl). Let u E IL^fl) and 
P(y, D)u E L 2 (fl). 

In fl we assume the existence of open connected subsets A*, flo,fc, a connected set T 
and functions ip k for k — 1, 2,..., K defined in this way: 

1. ip k E C 2,p (fl) for some p E (0,1), such that p(y , ip' k (y)) 0 and ip k (y ) 7 ^ 0 in fl 0 ,fc, 
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where p(y,£) = —£g + is the principal symbol of P. 

2 . Assume that: supp (u) D T = 0. Define Tb = fl 0 1 n T and for k > 2 set 
T* = no,m(U J<t A 3 -uT). 

Assume that there exist values ipmin,k < ipmax,k such that: ((supp (u) fl Qo,k) \ 
U<fe A i) c {y,Mv) < Y max,k } j and 0 {y G Q 0 }k ;ipk(y) max , k} C T fc . 

3. Assume that ip m in,k is such that A k = {y e (fi 0 ,fc \ T fc ); ip m in,k < ip k {y) < ip m ax,k } 
is nonempty, connected and satisfies dist(<9Q 0 ,fc, A*,) > 0. 

4. Assume that A = uf = 1 A fc is a connected set. 


Theorem 1.2. Under the conditions of Assumptions A2-A3, define the open set 
fli = (Jfc=i ,fc\T containing A. Then for every 0 < 9 < 1 we have 


IMIl 2 (A) < Ci61 

Moreover, for any m G (0,1] we get 

IMli/ 1 —(A) < C 161 






L'( Hi)' 


Mlrpcno 


('“(i + rfer))"'' 


i)' 


The constant ci 6 i is calculated in the proof. 


The dependency of the constant Ci 6 i from the geometric parameters of the prob¬ 
lems and from ip and 6 is described in Proposition 12.41 


2 Global Stability 

Notations. We start by introducing some notations and definitions used in the rest 
of the article: first we consider y = (t, x) G 1 x M n the time-space variable and call 
£ = (£o, 0 its Fourier dual variable. We remind that the exponential pseudodiffer¬ 
ential operator in Theorem 12.31 is defined as ! 2r v = J-fi o \ t e ~ e ^°^ 2T with 

J- and representing respectively the Fourier transform and its inverse. Then 
e~ e \ D °\ / 2r is an integral operator with kernel ,/2e . We consider a 

pseudo-differential operator A(D 0 ) with symbol a(£o) £ Cq°(M), 0 < a < 1 , sup¬ 
ported in |£ 0 | < 2 and equal to one in |£ 0 | < 1 . Hence we can write A(!3\Dq\/uS)v = 
J r { 0 Yt a (Y\^o\/^)J r t'^ 0 v and the integral kernel is ^ ). We will often 

work under the Assumption A 2 , where the symbol a is of Gevrey class. The smooth 
localizer h(y) is always supported in \y\ < 2 and equal to one in \y\ < 1 . 

The norm of the Sobolev space Hp is defined as ||w|| SjT = ||(|£| 2 + T 2 ) s / 2 J-y_^-u|| L 2 , and 
the space H s corresponds to the case r = 1. 

According to our notations the positive coefficients denoted by c x with x > 100 
are defined just once, independently on the variables / 1 , r, and they are calculated 
explicitly in terms of the coefficients of the operator (11,11) and the geometric parame¬ 
ters. This is essential to finally recover the value of C 160 and the radii R , r in Table l4~3l 
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A first step is the following lemma, proven in [6j, introducing a property often used 
in this section. 


Lemma 2.1. Let A(Dq) be a pseudo-differential operator with symbol a(£o), where 
a G Cq^R) is a smooth localizer supported in |£o| < 2 and equal one in |£ 0 | < 1. 

Assume that f(y) G C'o°(R", (^“(R*)), where 0 < a < 1. Then, for every p > 0, 

> 2, v G L 2 oc (R n+1 ) there are two constants c we , Ci 07 independent of p such that 

a) \\A((3iD 0 /p)f(y)(l - A(D 0 /p))v|| 0 < c 107 e~ Cl06fia ||v|| 0 . 

Moreover, if v G L 2 (R n+1 ) and h G C^°(R n+1 ) is a function such that h = 1 on 
supp(f), then 

b) \\A((3iD 0 /p)fhv\\o < \\f\\ oo \\A(D 0 /p)h(y)v\\ 0 + c w7 e~ ci0 ^ a ||/zu||o . 

If v G Hff c (WL n+1 ), m> 1, then the estimate above holds also in H m (M n+1 ), 

c ) \\A((3 1 D 0 /p)f(l - A(D 0 /p))v\\ m < cio8e _cl06MQ ||n|| m . 

Proof. See [ 6 ] for the entire proof. Here we remind how to obtain the coefficients. 
a) On the set supp[(l — p))a(fiifl/p)] one obtains |£q — £o|“ > (/i — 2p/f3f) a and 

the assumption /(f,.) G Gq (R t ) implies, uniformly in x on a compact set K C R n 
and for some C 3 = c 3 (a, K ), cn 7 = 0117 ( 0 ;, K) and C 106 = 0117(1 — 2//3i)"/4, 

|J-^ K i_ So) [/(t',a;)]| < C3 e- cll7 l«o-«°l Q < c 3 e~ 2ci (2.1) 

The coefficient c 3 = c 3 (a, K) is proportional to c±j, the Gevrey parameter of /, that 

is mm 

\D K (f(s))\ < c[f\ |k| + 1 ) |K|/ “, 5 G supp(/). 

We have c 3 = cijVol(supp (/)) and C 117 = l/(ec 3 ) a . We then estimate in the Fourier 
space the operator A((3iD 0 /p)f(-)(l—A(D 0 /p)), with c 10 7 = a( \ 1/a -^rcr)) 172 - 

- a(^o// i ))^eoM))llo < c ?o7 e_2ci06/i “lkllo • 


□ 

According to the splitting y = (t,x), the conormal bundle in R” +1 with respect to 
the foliation x =const is defined as: 

N*F := {(y,£) G T*R" +1 ; with £ = (&,£) and £0 = 0}. 

Its reduction to a subset K C R n+1 is Tk '■= {(?/,£) G T*K , £ 0 = 0}, 
and its fibre in y 0 is := {(yo,£) G N*F}. 

We then recall the concept of conormally strongly pseudoconvex function, alias strongly 
pseudoconvex function with respect to P 011 r yo ( [281 [29] ). 
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Definition 2.2. A C 2 real valued function 0 is conormally strongly pseudoconvex 
with respect to P at y 0 if 


Re{p, {p,0}}(j/ o ,O > 0 
on p(j/o,0 = °> 0^£GT yo ; 


( 2 . 2 ) 


{p(y,£ + irft(y)),p{y,£ + iT(t)'(y))}/{2iT) > 0 
for y = yo, 0^(6 T yo , such that p(y 0 , f + irfi'(y 0 )) = 0, r > 0. 


(2.3) 


In particular, for the wave operator (11.11) the conditions are void for non-characteristic 
surfaces 0 = const. As consequence one can state the following Theorem (Theorem 
2.1 in @). where the Carleman-type estimate by Tataru is recalled. 


Theorem 2.3. Let fl be an open subset o/Kx M n . Let P(y, D) be the wave operator 
(II.Ill , with gi k {x) E C 1 (fl) J h ?, q E C°(Ll). Let y 0 E Ll and if E C 2,p (fl) be real 


valued, for some fixed p E (0, 1), such that if'{yfi) 0 0 and S = {y;if(y) = 0} being 
an oriented hypersurface non-characteristic in y 0 . 

Consequently there is A > 1 such that <f>(y) = exp(A0) is a conormally strongly 
pseudoconvex function with respect to P at y 0 . 

Then there is a real valued quadratic polynomial f defined in (14. ip with proper a > 0, 
and a ball B^ 2 {yo) such that f(y) < fi(y) wheny E Bn 2 — {y 0 } and f(yo) = 4>{yo); and 
f being a conormally strongly pseudoconvex function with respect to P in Br 2 . This 
implies that there exist eo, To, c^t, C2,t, R, such that, for each small enough e < eo 
and large enough r > tq, we have 


e -e|D 0 | 2 /2r 



Here u E il/ oc (0), with P(y, D)u E L 2 (fil) and supp{u) C Bn(y o). 

This last estimate was used in [Bj to prove local stability of the unique continuation 
with explicit coefficients. We recall these results in the following proposition. 

Proposition 2.4. Let P the wave operator cnr Then, under the Assumption Al- 
A2, and using the result of Theorem \2.3\. there exist two positive radii R and r such 
that the local stability results (i.e. Lemma I '2. 61 and Theorems 1.2 in fid j) hold true in 
every point of 12 0 ; with the same parameters. 

Moreover, starting by the Assumptions Al, we are able to calculate all the constants 
involved in the local stability in a uniform way over The geometric parameters are 
constructed in Table (14.3p while the derived constants are in the proof of Theorems 
1.1 and 1.2 of Jfff. 

All the constants depend on : 

- the coefficients in (11.ip and their bounds: 


Is^IcqQob l^lc°(n 0 ), klc 0 (n 0 )) a h h, 


- the assumptions on the domains: 


dxst^n ^-1 {<9O 0 ,> 0, 
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- the non-characteristic condition and the non-vanishing condition upon if>': 


pi = min p(y,ij/), Q = min \ip'(y)\, 

y£flo 


- the norms of ip (see for notations) 

1^ Ic°(f2 0 )> 1^ lc 0 (O 0 )’ \W\max,n 0 i 

- the norms of the smooth localizers, in time-space and frequency, together with their 
Gevrey parameters: 


|Xl( s )|c r2 (R), |fr( s )lc 2 (R)> |®(£o)|c°(R)• 

We then need to reformulate Lemma 2.6 of [6j in the case of more general assumptions. 

Lemma 2.5. Under the assumption Al, let yo E Oo and ip be the quadratic polynomial 
ip(y) = f(y) — f(yo), with f defined in (14.11) . Let 0 < a < 1 and x(s) E Gq^“(M) be 
a localizer supported in [—85,(5] and equal 1 in [—75,5/2], Let p > 0, 5 > 0, be given 
constants, b E C^°(M n+1 ) and a E C“(M). Let A(D 0 ) be a pseudodifferential operator 
with symbol a. If 

U e h\b 2R ), Pu E L\B 2R ), \\A(D 0 /p)b((y - y 0 )/R)Pu\\ 0 < c A e~ » a , 

then for each r >0, there are constants Cn 0 , Ci 0 g such that 

|| e -^o| 2 / 2 r e r^ ) 6 ( ^o )p(j/jjD)u || o < max{c 110 , c A }e 2T5 ~ Cl09fJ,a max{ 1, \\Pu\\ L p B2R 

Using Lemma 12.51 we now reformulate Theorem 1.1. of [6] with more general 
assumptions. 

Lemma 2.6. Under the assumptions A1-A2, let y 0 E S = {y,ip(y) = 0} be an C 2,p - 
oriented hypersurface, which is non-characteristic in y 0 and with ip'(yo) 0. 

We also assume that u E H 1 (ft) is such that supp (u) D B 2R (y 0 ) C {y',ip(y) < 0}. 

Let b E G ( l / a (M n+1 ) be Gevrey functions with compact support, with 0 < a < 1. Then, 
for p > 1, if for some positive coefficients c B , cp, c A 

IMI hhb 2RIuq) ) < cu, \\Pu\\lpb 2R M) < c P , \\A(D 0 /p)b((y - y 0 )/R)Pu\\ 0 < c A e~ pa , 
then, there are constants C 150 , C 131 , C 132 independent of p, such that 

\\A(D 0 /u)b({y - y 0 )/r)u\\ H i < Ci 50 e _Cl32/i “ , Vw< p a /{3c 13 f). 


Moreover c± 3 i and ci 32 are independent of cp,cp,c A , while C 150 depends on them. 


The dependency of all the constants is as described in Proposition Wff 


Proof. The proof is identical to the one of Theorem 1.1. in [ 6 j. Th. 12.31 is used for 
the function x{t)^{^~ r ^) u that is supported in B R (y 0 )n{y, f>(y) < <f{yf)}U{y, —85 < 
<p(y ) < 5}. Moreover [P, x(<p)5(™)]u = [P, x(t))]K^-) u (since D(xb)u = bD(x)u), 
while in B R (y 0 ) one has x(t)^{^^) u = 














Here we have just to recalculate the related coefficients, distinguishing the ones de¬ 
pendent upon the parameters cu,cp , ca from the ones independent of them. 

We first list of coefficients independent of cp, cp, ca, but dependent on the Gevrey 
parameters of the localizers and from the geometric constants r, R, 5 (see Table (I4.3p h 


cix — ciw(l/«), C2X — l/(eA r cix)“, cioi — cioi(a), C102 — cio 2 (a, C101), 
Gig = 8c lx {a), B = S a c lx {a) : \d^{y)\ c o (BR) < c lls (R)(> 1), 

C114 = clMUXllUl + l^'| 4 C0/5 4 + W'\ 2 C»m 

C 115 = 4 ,A W\a> + l)(3 3 e- 3 /<5 3 )(l + lx'^o/ 5 2 ), 


, (1 — O') (i -°0 (1—«) 


27 TC 


Cl20 — ( \ td ) “ B “ lnS+1 , C121 — —T( 2 )(-j "C120, 


a In B 


a 


. 2 , 

a' 


C 122 — max(l, 4 cn 8 (i?)ci 21 , C\ X /R) > 1, C 123 ~ (eCi 22 )° < 1, C 128 — 7^^ c i23 < 1, 

C 110 = (ci22(8/3)T(l/a)/[ac^3(Q!Ci28) 1/(a_1) ]) i/2 , 
C 109 = min (a/ e 5 / 36 ,Ci 28 / 2 , 1 ) < 1, Ci 30 = 
ci 3 i = max{-\/2(16) 6 , (v / 2(16) 6 3^- 1 )V^)/ci 23 , ((16) 6 ^)/(3V2)} > 1, 

C135 = r a c 2 X ^ r , C 132 = min(ci 3 5 (r),ci 3 7 ), 

ClS7 = l CW2S ° ( ^fi Clxn 


Then the coefficients dependent on cp, Cp, Ca are: 


cue = 3max(c 2 T max{c n o,c A } 2 max{l, c P } 2 , cn 4 c^, cn 5 c^), 
eng = max(cn 6 , c^cn 2 (l + r 0 3 )(l + / 8 2 )e l25r °), 

Cl34 - C ^(( rClX )3 r ( a ) a ( r a C2X )l/a^ Q , Cl35 )l/(a-l)) 

C 136 — 2cioi\/cii3 f \J (s/ 5) 2 + 1 e Cl02S ^ds + 


+c 101 (^ 4 i±%ffi (2 

mm{l) c i 30 /2} y 


f+OO 


e~ y ' /2 dy'+ / e- cl02lxT dx') 

Jr 

C129 = max(ci 34 ,ci 36 ). 


We rename C 129 with ci 50 to underline its new dependencies. □ 

We now introduce the main assumptions to prove the global stability result. 

We recall that the support condition in Lemma 12.61 is not fulfilled everywhere for 
u e The idea is that at each step one applies the local stability result of 

Lemma 12.61 in a ball centred in the point y 3 and then one removes from supp (u) (in a 
smooth way) a part of the ball B r (yj) already calculated, for example by subtracting 
by 6 ( 2 ( 2 / — Vj)/ r ) u j, which is supported in B r (jjj). Then Uj + i fulfills the support 
condition in Lemma [2.61 in the ball B 2 R(yj. j-i), also due to our Assumption Al or A3. 
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Assumption A4 Let h2 0 and ip be as in Assumption Al. Then consider r and R 
the uniform radii defined in Proposition 12.41 

We define the set of points £ = {y 3 G O 0 ? j = 1, N}, such that Q a C UyLi BriUj) T 

Oo, in the following way: 

1. Let yi G h2o be the maximum point for ip in Q a . 

Set U\ — u and u 2 (y ) = (l — b( 2 ^ y / yi ^ ))ui, 

2. Let y 2 G O 0 be the maximum for ip in £l a \B r / 2 (yi), 

3. In general, let y 3 G fl 0 be the maximum for ip in Q a \ Ui=i bd r / 2 (ykP), he. y 3 G 

argma x{ip(y), y G (fi a \ Lid B r/2 (y k ))}. 

Then we dehne: 

u j = JJ (1 -b k )u, b k \=b(^-—^-). (2.4) 

k=l..j-l 

Each y 3 lies on the surface Sj = {y, ip(y) = ip^yf)}. Notice that, since | y 3 —y k \ > r /2 
for j 0 k, 


N < Ci70 


Voi(n 0 


^n +1 ( 


AyJ max(6i ,1) 


n+1 ! 


(2.5) 


where oj n+ 1 is the volume of the ball of radius one in R n+1 , where we consider the 
following bound for the coefficients a\5^ k < g jk (x) < b\S jk . 

We finally dehne l(y) G Qj’°(M n+1 ) a localizer such that l — 1 on C — {y G 
h2 0 ; dist(?/, <9f2o) > 0 < l < 1 and supp(Z) C Oq• Observe that U k=1 B 2 n(y k ) C £. 


We now can formulate a stability estimate of inverse exponential type for the low 
temporal frequencies of Uj. 

Theorem 2.7. Under the Assumptions A1-A2-A/, lety k G £ and let b G 
be a Gevrey functions of class 1/a with compact support, such that 0 < a < 1. 

Then, there exist constants R, r with R > 2r > 0, and Ci 59 > 1 such that for y > ci 59 
there are coefficients C 151 , c 152 , c 154 , c 455 , Ci 56 , /3, N for which, 

if 

IM|/p(nO = 1, \\Pu\\ L 2 {ni) < 1, \\A(^)l(y)Pu \\ L 2 < exp (-//*), (2.6) 

then calling Hi — y and y 3 = Ci 56 /r“__ 1 for 2 < j < N , we have fij > 1 and 

\\ u j\\m(B 2R ( yj )) < Ci 52 , \\Puj\\L2(B 2R ( yj )) < C 153 , (2.7) 

\\A(D 0 /(ijWiy - y j )/R)Pu j \\ 0 < c 15iJ exp(-/i“), (2.8) 

and consequently 

\\A(D 0 /uj)b((y - y^/r'jUjWui < c 155tj exp(-c 132 yf), Vo; < y? / (3c 131 ). (2.9) 

The radii r and R are defined in Table (14.31) . while the coefficients c k are calculated 
in the proof of the Theorem. 
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Proof. Let b G Gg^“(M n+1 ) be a localizer with support as in Assumption 2. 

Observe that according to onr definitions we have : 

Briyj) C supp b((y - y 3 )/r) C B 2r (y 3 ) C B R {y 3 ) C supp b((y - %)/A) C B 2R (yj ), 

We now proceed step by step. 

Step 1. We consider y\ G £ defined in Assumption A4. From the hypotheses (12.61) 
the following inequalities hold true for u\ = u: 

IMI^OWs/i)) < 1, \\Pu\\L2(B 2R ( yi )) < 1- 

From the definition of l in Assumption A4 and applying Lemma 12.11 (b) with / = 
h = l, /3i = (3, y = /3y we get 

||A(—)fe( y yi )%)Pu || L 2 < ||A(^)/(?/)Pu || L 2 +c W7 exp(-c 106 /3 a y a )\\l(y)Pu \\ L 2 

y K py 

< exp(—//“) + Ci 07 exp(-^“) < Ci 54 i iexp(-/i“) (2.10) 
with C 1541 = 1 + C 107 and where j3 > 2 is a parameter chosen as: 

/3 = 2 + (J^) V “ (2.11) 

Cil7 

in order to have Cioe/3 a = 1. Indeed, applying Lemma 12.11 (b). one gets C 106 = 
^(1 - |)“ where c 117 = 1/(2c 3 )“ and c 3 = ci, b (P) • Vol(supp (b(^jp-))), with c hb (R) 
the Gevrey parameter associated with 6 (|;). For the calculation of C 107 see lemma 
O Notice that C 106 and C 107 are independent of y \, since the calculation is invariant 
up to translations. 

Calling ip(y) = if(y) — ip(yi) we notice that u fulfils supp (u) fl B 2R (y) C {y',if(y) < 
if{yi)} We are then allowed to apply Lemma [2761 with y 0 = yi, if = -0, cu = 1, Cp = 

1 ,ca = Ci 54 ; i and calling <+ 55 ^ = C 150 , 

11 a / D» ./y 2/1 \ n / / \j ^ y 

A(—) 6 (- )u\\ H 1 < C 155 1 exp(— 0132/1 ), Vca < --. 

a; r 3 ci 3 i 

Step j > 1. 

Here we consider y 3 G £ and u 3 defined in (I2.4[) and notice that supp {u 3 ) C supp(w)\U^~ 1 
B r / 2 (y k ) and that u 3 = u on supp(u)\ B r {y k ). 

Calling ip(y) = ip(y) — tf(yj) we notice that by construction Uj is such that supp (u 3 )C\ 
B 2R (y 3 ) C {y',if(y) < ip(yj)}. We then will apply Lemma [2761 with ip — ip and 
yo = Vj- 

We start by calculating the first estimate in (12.71) : 

\\ u j\\H^{B 2R { yj )) < \\u\\m(B 2R { yj )) + |V (l -b( ^ ^ yk ^))\ c o\\u\\ L 2( B2R{y .){2.12) 

k=0..j—l 

< 2(1 + j\—^-)\\u\\ H i(B 2R ( yj ))- 
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Since j < N we get a uniform bound for all j 


\V\ 


i ^ i C® 

\\ U j\\H 1 (B 2R (y j )) < Cl52, C 152 = 2(1 + N —— ). 
Then we consider the second estimate in (12.71) 


(2.13) 


\\ Pu j\\L 2 (B 2R ( yj )) < 11 Pu\\ l 2 (B2R(yj )) + II [-P, ]^[ (l - —))} u \\L^B 2R (yj^-^) 


< 1 + 2j(l + n 2 \g kr \ c o + \h s \ c o)(^l + 


k=0..j—l 

W\c° , \b"\c° 


+ (j - 1) 


I* 


/12 
I CO' 


rj^2 


\ U \\h' (B 2Jl ( W )) - C 153, 


where the commutator is, for b k = b(2(y — y k )/r ): 

[6fc, -P]« = ( -P 2 bk)u + 2D 0 b k D 0 u - 2g hr {x)D Xh b k D Xr u + ih s {x)D Xs b k u 
with P 2 = —-Dq + 9 hr ( x )Dh.D r and, for all j < AC, 


C 153 — 1 + 2N(1 + n 2 |5' fcr |c'0 + I /i S | O ' 0 ) (- — - h 


Ife'lco , \b"\c° 


+ (i V-l)E|i), (2,15) 


The third estimate (I2.8[) requires information of Step j — 1. 

Like in (12. lOj) . from the definition of l and applying Lemma \2 .11 (b) with / = 
h = l, f3i = (3, y, = PjAj we get 


R yj )Ky) Pu j\\L* < \\A{^)l(y)P Uj \\ L , + Ci 5 3 Cio 7 exp(-/i"). (2.16) 

where j3 > 2 is the parameter (12.111) . 

The first term on the right hand side of (12.161) becomes 

WM^Viv) Pu ib < \\A(9±)l(y)Pu ,^||„ + IW^Wi-i^ViHo 

+W A ^)[bi-uP]u i -i\W (2.17) 

One can find recursively the estimate above for j = 1 by using (12. 6 j) with ci 62 ,i = 1, 
and stating for j — 1 


\\M 


Do 

ft 9 j —i 


■)l(y)Puj- i||o < ci 62 j-ie 0“-i, 


(2.18) 


with Ci 62 j_i a positive parameter. 

By the inductive hypothesis and in analogy with (I2.10p . 

II A(-^-)b(- — j^—)l{y)Puj-i\\o < ci54j-ie“^-i, 

l^j — 1 JA' 

where c 154 ,j-i = ci 62 j-i + c 153 c 107 . 

The first term on the right hand side of (12.171) becomes, for gj < [ij_ i/2, 

P(# L )f(y)*S-illo < illo- 

PPj PPj-i 


(2.19) 
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The second term on the right hand side of (12.171) becomes, for 2/n, < 

\\A(^-)l(y)b j -lPu j - 1 \\ 0 < || M q D ° )Kv) Pu j- l||o + Cl53Cl64exp(-Ci 65 /i“_i) 

PUj P/ij-i 

where by lemma [2711 b) with / = bj- 1 , h = l, f3 x = 3, n = /3fJ,j- 1 , we have Ci 64 = C 107 , 
Ci 65 = Cio 6 /d“ = Cn7/3 a /(3 Q 4). Notice that Ci 65 and Ci 64 are independent of y 3 , since 
the calculation is invariant up to translations. 

The term with the commutator in (12. 1 7p can be split in the following way: 

I 1 + I 2 + I 3 = ||7l( ——) ((— Pibj-i) + ih s (x)D Xs bj-i)Uj-i\\L 2 

PUj 

+ ||-^(o—) { 2 D 0 bj-iD 0 Uj-i )\\ L 2 
Pflj 

+\\A(-^-)( 2 g k '-(x)D n b t . 1 D Xr u j . 1 )\\ L ,. 

We notice that the localizer b (—— — 1 - ) — 1 on t .P y. then we multiply 

u,j _i in /] with it to keep its support in B 2 r (yj_ 1 ) in order to use the estimates of 
H°-_ x 

Step j — 1. For z/ < —— a positive parameter, one has 

3Ci3i 

+\\A(^-)((-P 2 b j p+ih’(x)D„b j . 1 )(l - 

2 

< ci 55j -i| - P2&J-1 + ^(a:) J Dx s ^-i|c o exp(-ci 32 /i“_ 1 ) 
TC107C152 (1 + n 2 |5f fcr | c o + |/r s | c o)exp(-ci 06 ^“) 

Notice that the first estimate on the right hand side is done by using the inductive 
hypothesis and by applying to the term \\A(Bj)b p J — ! ! )u./-i||u 2 Lemma 12.61 with 
coefficients cp = Ci 52 , cp = c 153 , ca = Ci 54j _i defined in (12.19j) and then calling 
Ci 55 j-i the resulting coefficient Ci 50 = ci 50 (ci 52 , Ci 53 , Cimj-i). 

For the second term on the right hand side we assume that 2 /3/j.j < u/3 in order to 
write, both with s = 0 (i.e. L 2 ) and s = 1 (i.e H 1 ): 

< Il4(~~)u||s- (2.20) 

PUj v 

Then we apply Lemma [2.11 (a) with j3 x = 3, n = v and / of this form (after moving 
out of the norm g kh , h s and the complex variable) 

fi(t, x) = d%bj _i + d Xr d Xh bj -1 + d Xt bj- 1 , (2.21) 

involving just derivatives of smooth functions in Co°(K", Gc/ cl (Mt)). 

To recover an expressions for the coefficients we recall that the k 2 —derivative of 
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h E G X J a (with Gevrey constant Ci^) is: 

\D K \D^h{s))\ < t >i | + |,e2| + 1 (| /Cl | + \ K2 \ + 1 )(l"il+M)/« 


( 2 . 22 ) 


< c ^l + l' s il+ 1 2 l« 2 l(l«il+l"2|)/a e |i*2||*i|/aQ Ki | + 1 )|*i|/« j s G supp (ft). 

In our case we must just consider time derivatives, in order to estimate (12.11) . Since 
the translations play no role for the Fourier transform, we can calculate coefficients 
independently upon j. Call c P 6, Ci^, c\pi the Gevrey coefficients of the functions 
b(y), D x b(y ), D 2 x b(y). Then, define 

c h = 2 2/a e 2/ “c 2 fe (r) + a , 6 »(r) + c p6 /(r). 

Analogously we can get the values for the functions associated with A, J 2 (see below 
for definition of / 2 , /a): 

c/ 2 = 2 1/a+1 e 1/a ci i b(r) + 2 2/ “ +1 e 2/a c p& (r), c h = 4ci >6 /(r) + 2ci,y/(r). 

In analogy to the computations above we can calculate cd/ 2 , cd/ 3 (the Gevrey param¬ 
eters of Df 2 = d t (2d t bj-i) + V x (2d t b j _ 1 ) , L>/ 3 = d t {2d Xk b 3 _i)) + V x (2^6 i _i)), in 
order to apply Lemma [2.11 (c) with H 1 — norms. 

Now call c comm = Cf x +Cf 2 + Cf 3 + CDf 2 + CDf 3 the biggest Gevrey parameter, common to 
all the functions /j, / 2 , / 3 inside the commutator, set c 3 = c comm ■ maxj Vol(supp (/*)), 
and Cn 7 = l/(ec 3 )". Then, define the following coefficients in Lemma 12.11 that are 
independent of the center point yy. 


1 , 1 1 

c i06 ~ TTFnT-w Cl07 — Cl08 — f—J —7--T 

(3°4)(ec 3 )“ 3 y a'a(cm) 1/a (ac 106 ) 


a/2 


1 

a —1 


(2.23) 


Next we estimate J 2 moving the derivative D 0 of u 3 _\ in front of the integrand, 
then multiplying u 3 _\ with b( y >h ~ 1 ), and finally adding and subtracting operators 

A(D 0 /u) with v < 

jt/ jU/j l' / 

+ ||A(^) J D 0 (2 J Do^-i) [A(^) + (1 - A(^))]6( (y ~^ l) )]u J _ 1 |l L2 




v 


2 

< c 155j-l 12ZT 0 />J — 1 \ci exp(—Cl 32 /ij_ 1 ) + Ci5 2 Ci08 exp( —Ci 06 ^“) 

2 

+ci55,j-i|-Do(2_Dofrj-i)|c , ° exp(—ci 32 /r“_!) + cis 2 cio7 exp(— cia§v a ) 

To get the estimate above we apply twice Lemma [2761 with the same parameters as in 
I\. Next using (I2.20j) we estimate the terms \\A(^ 3 -)f(l — A( — ))u|| s using Lemma 
12.11 c) and a). 

Proceeding like with ii, we have to calculate the time-Fourier transform of: 


h{y) = Zdtbj-i + 2 d 2 b J _ ll 
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and the associated coefficients are (12.231) . Finally, moving the derivative D Xr of Uj-\ 
in front of the integrand, multiplying with b( u / °~ 1 ), adding and subtracting 

operators A(D 0 /u) with u < 

h < \\A(^-)2g kr (x)D xl ,b i . 1 [A(^-) + (1 - A A)]6( (y ~ 

|L/ f-v j L' I 

+\\A(^-)D Ir (2g kr (x)D xt b j _ 1 )[AA + (1 - A(^))]b( ^ |U, 

fS j IS IS / 

< C155j-i 1 2 ng kr D k bj~i | c 1 exp(-ci 32 /i“l 1 ) + c lb 2 c W sn 2 \g kr \ c ^ exp(-ci 06 i / “) 
+c 1 55,j-i\D r {2g kr D k b j _ 1 )\co exp(-ci 32 /i“! 1 ) + c 1 Q 7 c lb 2 n 2 \g kr \ c ^ exp(-c 106 u a ). 

Proceeding like with /, we have to calculate the time-Fourier transform of 
Mv) = 29x fc 6j-i + 2dx k bj~i + 2dl k bj-i. 

and the associated coefficients are (I2.23p . By collecting all the terms of the estimate 
for (12.171) . the bound for (12.161) becomes 

11 A(—)b( V y3 )Puj\\ 0 < \\A(^-)l(y)Puj \\ L 2 + ci 53 cio 7 exp(-/i“) 

/ij K ppj 

< ci 62 j (max (exp(-/i"_ 1 ), exp(-ci 65/ n“_ 1 ), exp(-ci 32 /x“l 1 ), exp(-ci 06 i/“))) 

+Ci5 3 Cio7exp(— fij) (2.24) 

where, for all j > 2, 

c 162 j = 2ci6 2 y-i + Ci5 3 Ci64 + Ci55j_i| — P 2 bj_i + h s (x)D Xs bj _\| c>o 

+C107C152(1 + n 2 \g kr \ c o + \h s \ c o) 

"PCl55j_i \2D(jbj—i \qi + Ci5 2 Ci08 
+ c i55j-i I D 0 (2D 0 bj_i) | C° + C152C107 

+Ci55,j-i\2ng kr Dkbj-i\ci + ci5 2 cio8^ 2 |fi ,Ar |c 1 

+C 155 ,j-i\D r (2g kr Dkbj-i)\c° + cio7Cis 2 n g kr \c 1 - 

In order to write (j2.24[) in the form 

II A(—)b( V yj )Pttjl| 0 < ci 54j e _ ^, 

jfl j JlL 

we set in (12.8[) 


C 154J — C 162J + Ci5 3 Ci07 


and we look for pj of the form pj = Ci 56j such that, for /3 as in (12.lip , and 
collecting all the constraints on pj used in the proof, 




c 156 — c 156,i — mm 


1 tf- 


4-1 


6/3 3ci 3 i 
1 

18/3ci 3 i ’ 


and consider 

1/a 

_V a l/“ C 106 

-132 j c 165 1 


3c 


131 


O' 


(2.25) 


15 




























The right hand side of (12.25ft is independent upon j due to the definition of c 165 , Ci 0e 
in (12.231) , and the fact that Ci 32 and Ci 3 i do not change during the iteration (see proof 
of Lemma [2.61) . Therefore we define a parameter Ci 5 g(< 1) independent of j by the 
formula (j2.25j) . We then estimate /ij from below 

/L = c 156///-l = Ci56(ci56/i“_ 2 ) Q = <H56(uL56(ci56/d/_ 3 )“)“ 

_ l+a+a 2 +...+ch -2 a J_1 \ 1/(1—“) ,,ch _1 

— c 156 r '156 r 

where we apply 1 < Y7m=o ^ — 1/(1 — o>). Finally, in order to obtain the requested 
condition (ij > 1 , we set 

_ i 

c i 56 1 ^°V qJ 1 > 1; Vj G [l,iV], which implies /i > c 15 

i 

and we finally find C 159 = c 15 g JV_1(1_ “ ) in the line before (12.61) . 

By applying Lemma ESI with c v = Ci 52 , c P = Ci 53 , c A = c 154 j, 

c i 55 ,i = c i 5 o( c i 52 , C 153 , Ci 54 j), one obtains the last inequality (12.9[) of the Theorem. □ 


Remark 2.8. 1. In order to work with the pseudodifferential operators e.~ e \ D ^ P T 
and A(/3\Do\/u;) one needs smooth functions in the time variable. Hence one should 
first operate a proper regularization in the time variable. We proceed in the same 
way as done in [10] or H3. Observe that the functions u, Pu and Uj, Puj are always 
multiplied by a smooth localizer when A(Dq) and e~ D ° are applied to them. 

2. About the construction. 

a) Assumption Al (and analogously A3) implies that: (supp {u) fl L2 0 ) C {y^iv) < 
0 max }; a nd that the level sets {y G 9 O ;00/) = c}, with c G [0 m m,0ma0, are con¬ 
tained in T U £l a . An example of this construction is in section [3] 

b) Assumption Al and A3 can be relaxed in this way. Instead of defining fl a (or 
Afc), we just observe that the assumptions on 0, f2 0 together with (supp (u ) 0 flo) C 
{r/;0(r/) < VWr} imply the existence of a non empty set fl a C f2 0 f° r which Theo¬ 
rem [TTT] holds. can be defined as U jB r (yj), with y g G £ (see Assumption 4) such 
that the support condition supp (uj) fl B 2 R(yj ) C {y G Ho) 0 < 0(?/j)} is fulfilled 
for every j. This construction requires to follow step by step the local iteration and 
sometimes this is difficult. That is why the a-priori knowledge that the level set 
{y G f2 0 \ T;0(?/) = 0 min } is strictly contained in f2 0 is useful, even if it excludes for 
example the case where the level sets of 0 are parallel hyperplanes and supp (w) is 
on one side of one of them. 

3. In Theorem 12.71 we have worked under the assumptions 

IMIrrqni) = 1, \\Pu\\ L 2 {ni) < 1, \\A(D 0 /(Pfj))l(y)Pu \\ L 2 < exp(-//“), 

in order to apply Theorem 11.11 easily. One can generalize the assumptions by setting 
IMItfi(ni) < cu, g , ||Pu||i, 2 (ni) < Cpg, \\A(D Q /(/3n))l(y)Pu\\ L 2 < c Ag exp (-//*), 
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and by changing the coefficients Ci 52 , Ci 53 , Ci 54 j, Ci 55 j accordingly. 

This gives a statement of global stability of the unique continuation for low temporal 
frequencies. 

4. Notice that Lemma [2751 Lemma [2761 and Theorems 12.41 [Lilli.21 can be reformulated 
for localizers supported on cylinders (instead of on balls), defined on cylinders 

C s (y 0 ) = {(t, x) : \t - t 0 1 < s, \x~x 0 \< s}, 


by observing that: 


C r /V 2 (yo) C B r (yo) C B R (y 0 ) C C^ R (y 0 ). 

The advantage is to be able to reduce the assumptions on the regularity of the 
x— localizers. Namely, one can replace b(y) in G ( l / Q (R n+1 ) with the product b fl (t)b sp (x), 
where b u G G,!/“(]R t ) and b sp G Gq(IRC) are supported in £> 2 , equal to one in B\ and 
0 ^ bu , b S p ^ 1. 

One can also replace the global localizer l in the proof of Theorems 12.41 with l(y) = 
J2m=i ki,m(t)lsp,m(x), with l = 1 in {y G O 0 ; dist(y,dQ 0 ) > , which contains 

u k=i C V 2 R(yk)^ and supp(Z) C O 0 , and l ti>m G G^(M) and l sPt m G Gp(M n ). 

2.1 Proof of Theorem 11.11 

Proof. We consider two cases: 

1 1 

Case A. Assume ||Pu||z, 2 (ni ) > exp(—ci 59 ) ||rt|| (f2l) , where c 159 := c 15 ( 6 1_Q ° “ JV_1 > 1 
has been defined in (12.251) . Then the estimate is trivial 

IMMn.) < ll«llffi(Oi) < (ln(l + exp(ci 59 )))*-- ~ e ■ 

o i+ fe) 

Case B. Assume now H-P^Hl 2 ^!) < exp(—Ci 5 g)||tt||iji(n 1 ) and without restriction of 
generality take IMIn^Oi) = 1. Our aim is to consider separately estimates for low and 
high temporal frequencies. Let A (Do) be a pseudo-differential operator with symbol 
a(£ 0 ) G Gp /a , defined in Assumption A2. Let b G G, l / Q (R n+1 ) be another localizer 
with support like in Assumption A2. 

The parameter a G (0,1) is then common to all the localizers in time, space and 
temporal frequency. Let yj G £ be the set of the points defined in Assumption A4 
and consider the balls B r (y 3 ) centred in those points. Observe that according to our 
definitions we have : 

B r (Vj) C supp b((y - yf)/r) C B 2r (vj) C B R (yj) C supp b((y - yj)/R ) C B 2R {y j ), 

Recall that b(^-^-) = 1 in B r (yj) and observe that u 3 = u in B r (y.,) \U J S I [ £> r .(y. 5 ), with 
Uj defined in (12.41) . 

We then cover Q a by the disjoint sets B r (yj)\ulz\B r (y s ) and operate the initial 
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estimate as follows: 


M IU 2 (n«) ^ \\ U \\^{B r { yi )) + \\ U \\^(B r (y 2 )\B r ( yi )) + W U \\ L ^B r (y 3 )\U 2 s=1 B r (y s )) 

+ ••• + H U llL 2 (i? r (j/ i v)\uf = - 1 1 i?.(^))( 2 - 26 ) 

= ||6(^ r Vl )u\\ L 2(B r (y i)) + \\b( J r V ~)u2\\L 2 {B r (y2)\B r ( yi )) 

+ ... + || b( V r JN )u N \\ L2(Br {yN )\ u ^ Br (y 3 )) 


N n N 

< E M(—W ?L - + E IK 1 - A (— ))K ILJ!i )«jllw := Hi + Hi. 

oj r z —■' ' u r 


In the last estimate we have chosen cu > 0 and split all the terms in their low and 
high temporal component, i.e. 

< \\A(^-)b( y -^i ) Ui ||y + ||(i - a(^))6(^)%|U.. 

To estimate IL> we have 

||(1 - < ||(1 - a(f„/^))-WK(?/ - y^Mufymh (2.27) 

< A f [ - y^Mujfax^dxdto < \\\b((y - Vj)/r)uj{y)f Hl . 

u J |£ 0 |>w>/R n w 

To estimate H i we first consider /i > C 159 and we set \\Pu\\l 2 ^^ = e _/i , that implies 
\\A(^)l(y)Pu\\o < e _M “, for all ( > 0. Then we choose ui = n%/(3ci3i) and apply 
Theorem 12.71 to each term of the sum: 

\\A(^-)b( y r Vj )uj\\ L 2 < ci 55i Arexp(-ci32Ai^ 2 ). (2.28) 

This is possible since > 1 is the smallest time frequency of the set fij, while Ci 55) tv 
is the biggest coefficient Ci 55 j, j = l,..,N. 

Collecting the two bounds and reminding that Hn > c i 56 V" 1 > 1 , we have: 


3c N 

M|l/*(n„) < Nc 155 :N exp( Ci 32 /i^) + - -^Ew 

,lN j =1 

< Nc 155>n exp (- C132 (c}/i QjV_1 )" 2 ) 
3iVci3i ( l^lc 0 ^ 

,l/(l-c) Q iv —1 so \ + —:— ) C 152 


y-yj 


)% I Iff 1 (Hi) 


< 


( C 156 '^ 

c 158 c 158 


< 


2“ ci 58 ||'u||.fl-i(n 1 ) 

^ f—■ln(||PM||x, 2 ( n 1 )))“ iV " (In (1 + )V 

V V \\ Pu \\ L 2(n 1 ) J J 

where C 156 is defined in (12.25ft and 

Cl58 = Nci55 t N + 31VC 131 C 152 (l + -—— )c 15 / < ' ^. 
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In the last step we have applied In (y) > ln(l + y)/2 for y = ||M||_H'i(n 1 )/||-Pu||L 2 (Qi) > e, 
and then we have returned to the original notation. 

Now we choose a such that a = (9) l t N and which belongs to (0,1) so that, defining 
ci 60 = (ln(l + e Cl59(e) )) e + 2% 58 (6 i ), we obtain the result. □ 

In the previous theorem the dependency of Ci 60 upon 6 is very bad. 

For some applications it is better to keep a and N independent and formulate the 
following consequence: 

Corollary 2.9. Consider the assumptions of Theorem, li.il Then for every 0 < a < 1 
we have 


M|L2(Q a ) < Ci 60 


\u 




(In (l + 


u ffi 


g‘(n i) 


■))' 


11 ^ 111 , 2 ( 0 !)' 

with ci 6 o = (ln(l + e ci59 ))“ + 2 qJV ci 58 . Here N < c 170 given by (12.5ft . 


2.2 Proof of Theorem 11.2 



FIGURE 1. A possible construction of the domains T,Tlo tk , A k 

Proof. Initialization of the radii and the localizers: Let A{Df) be a pseudo-differential 
operator with symbol a(£o) £ Cq/“(M), dehned in Assumption A2. We define the 
localizer b(y) e (K n+1 ) with support like in Assumption A2. 

Using Assumption A3, in each domain Oo,fc we can calculate a table like (14. 3 j) . with 
f2 0 ,k i n place of O 0 , and A^ in place of f l a , and where all the constants dependency 
is described in Proposition 12.41 By comparing the tables of the several we can 

_ j2 

consider R 2 = mm k R 2 ^ and find R = ^(16 + R 2 the common radius for the 

local stability in A*. After fixing R , we reduce also the values r k so that r = min k r k 
is the common radius of the ball where the L 2 local estimate can be performed in A fc . 
Construction of the set E and the functions uf For f2 0 = ^o,i) = A] and = ifi, 

we define y 3 6 £\ the set of the maximal points for -i/q, according to the procedure 
in Assumption A4. Call N\ the number of points of the covering of Ai, i.e. Ai C 


19 














Uj^iThen we remove Ax from £l and we restart the procedure with the set 

A 2 - 

Namely, for h2 0 = O 0 , 2 \Ai, V’ = "02 and = A 2 , we define G £ 2 according to the 
procedure of Assumption A4, where we use the indexing j = N\ + 1, N 2 . 

When also A 2 is covered, one skips to A 3 and so on. 

At the end one can define the set of points £ = Lif =1 £ k of cardinality N = 1 -^4- 

Consider the balls ) centred in those points and define Uj as in (12.4f . 

Observe that according to our definitions we have : 

B r (Vj) c supp b((y - yj)/r) C B 2r (yj) C B R {yj) C supp b((y - yj)/R ) C B 2R {yj ), 

Moreover 6(^ Z ^ L ) = 1 in B r (yj) and observe that Uj = u in Z3 r () \ I j Z? r (y 5 ), with 
Uj defined in (12.4|) . 

Uniform parameters: With the previous assumptions Lemma 12.61 can then be applied 
in each ball B 2R (yj ), with yj G £ k in place of y 0 , and we call Ci 3 i )fe , 032,fc, 065,fc, Oso,fc 
the related parameters, that are constant for every ball in the region O 0 ,fc: but in 
principle change from region to region. 

Therefore the following uniform constants are introduced: 


031 ,* 


max C 131 k , O32,* — m in 032,fc, 065,* 

k=l..K k=l..K 


min 

k=l..K 


065 ,ki 


cise,* = mm 


18/3o3i,* 


1/a 1/a 

032 ,*) 065 ,*; 


T/ Q 

"106 


3c 


131,* 




(2.29) 


We define l(y) G Cq°(M" +1 ) a localizer such that 0 < l < 1, l — 1 on the set 
{y e Ufc=i^o,fc; dist(y,dQ, 0 ) > x); which contains U f =1 B 2R (yj), and supp(Z) C fi 0 . 
In particular we consider /? as in ( 12 . lip , and the related Ci 0 6 - 
Construction: 

We consider two cases: 

1 1 

Case A. Assume ||Pu|| L 2 (Ql) > exp(-ci 59 i *)||tt||Hi ( n 1 ), where Ci 59j * := c 15 ff a) “ JV_1 > 1. 
Then the estimate is trivial 


IMU 2 ( A) < IMMo) < (ln(l + exp(ci 5 9 ))) 0 - 

w+fer 

Case B. Assume now H-Pwllzp^) < exp(— 059 ,*)||w||jj-i(n 1 ) and without restriction of 
generality take |M|,fp(ni) = 1 - Our aim is to consider separately estimates for low 
and high temporal frequencies. We cover A by the disjoint sets B r (yj) \Lk s I \B r (y s ) 
and operate the initial estimate as follows: 


I MU 2 (A) < ||w||L2(B r .(i/i)) + 11U11L 2 (B r (y 2 )\B r (j/i)) + ••• + ll u lll, 2 (B r (j)iv)\U^ = q 1 S r ({/ s )) 

= l^( ~ )' U IU 2 (0(j/i)) + ••• + \\K ~ ) u N\\i?(B r (j/ivlXu^Ti 1 B r (y s )) 

N n - N n - 

< E \\ A (£-)i>C-f^)uj\\Li + E1K 1 - := + h 2 . 

j =1 3 = 1 


20 


















In the last estimate we took oj > 0 and split all the terms in their low and high 
temporal component, i.e. 

+ ll(i - aQ)) 6(^) Uj ||y. 

To estimate H 2 we have 

||(1 -A{D Q /u))b{{y-y j )/r)u j {y)f I? < ^WKiv ~ Vj)/ r )uj{y)\\m- 

To estimate H 1 we first observe that supp(6(^-^ L )Mj) is in f2 0 ,i for j — 1,Ah, where 
ipx is dehned. Then supp (b(^-^-)uj) is in flog for j = Ai + 1, N 2 where ^2 is defined, 
and so on. 

Consider /i > ci 59j * and we set ||Pw||z 2 ( 0l ) = e _M , that implies \\A(^)l(y)Pu\\ 0 < 
e _At “, for all ( > 0 . 

We also set c 152 = 2(1 + 1V^), c 153 = 1 + 2A(1 +n 2 \g kr \ c o + \h s \ c o)(^ + ^ + 

\b'\ 2 

(A — 1) r 2 ° ), and )ij = c 156 ) */i“_,. Then we choose lu = /( 30131 ,*) and apply 

Theorem 12.71 to each term of the sum: 

, r ^ IU 2 < Ci5 5i jvexp(-ci32,*^ 2 ). 

This is possible since /tat > 1 is the smallest time frequency of the set fij, while 0155 , at 
is the biggest coefficient C 155 J, j = l,..,N. 

Collecting the two bounds and reminding that /j,n > “V ^ 1 > 1 , we have: 

Q N 

II II / ( a 2 \ , oCi3i >!t ^ .y — ]Jj. 

|M|i 2 (A) < tVCi5 5i Arexp(— Ci32,*I^N ) H-q— / y WK - )Uj || H 1 ^) 

/% j=1 r 

na N || || 

< Cl58,* _ _ Cl58,* _ < 1 Cl58,* 11^ 11 1 (Oi) 

H aN (-ln(||Pu|| L 2 (ni )))^ ~ ( ln ( 1+ )) aAr ’ 

V V H P “lh2 (ni) // 

where 

c 158,* — ^Ci55,AT + 3ATci3i j<t Ci52 (l + | ) C 156 ^ \ 

Now we choose a such that a = {6) l ^ N and which belongs to (0,1) so that, defining 
ci 61 = (ln(l + e Cl59 -* (0) )) e + 2 e Ci 58 ,*( 0 )) we obtain the result. □ 


3 Applications 

Assumption A5 Assume that M = M n and on M we have a metric tensor g 
satisfying 

oq! < [9jk{x)] n j yk =i < b 0 1, and \\gjk\\c*{M) < h- (3.1) 
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Note that then 


b^I < ^‘(i)]" t=1 < a 0 -n, (3.2) 

||5 ,jfc ||c 1 (M) < h^o 2 ) 

y/a^d R n(x 1 ,x 2 ) < d g (x i,x 2 ) < •\/MRn(xi,^2), 

^llfllr < II^IU < -Lr|| r „, ( 6 i;*M, 

\/ao|Mk" < llhlls < tAolMk"; h e T X M. 

Here we assume that ao < 1 < bo- In Appendix A we call a\ = 1 and 61 = Og 1 - 

Note that (18.ip implies that 

I Sec | < A M = A M (a 0 , b 0 , b 3 ), 

where Sec is the sectional curvature of (M, g ). Also assume that the injectivity radius 

1 /o. 

of ( M,g ) satisfies inj (M,g) > i 0 with 0 < io < min(l, 7 r/A^ ). 

Consider the wave operator ( 11 . 11 ) . Assume that the lower order coefficients are such 
that 

lk J l|c 0 (M) + IMIc°(m) < b 3 . 

We fix the three positive parameters £, T, 7 as follows: 

i < i 0 /4, T >£, 0 < 7 < T — C 


In this section, we use the following dehnitions. 

Definition 3 . 1 . B g (z,ri) C M = M n is the ball with center z and radius r 1, defined 
using the Riemannian metric g. Also, B^n(x, ri) is the Euclidean ball in M n . For 
y — (t, x) G M x M, let 

C g (y, ri) = (t-r 1 ,t + n) x B g (x, 77 ) ( 3 . 3 ) 

andC(|/, ri) = C^n+i(y, ri) = (t — 77, £ + 77) x £>Rn(x,ri). Also, denote 
dRx(wr-, g )((ti,xi), (t 2 ,x 2 )) = max(|H - t 2 \,d g (x 1 ,x 2 )) 

and 


diRx(K-,e)((n,xi), (t 2 ,x 2 )) = max(|R - t 2 \,dmn(xi,x 2 )). 

Let z G M, and define the hyperbolic function as and 

-0(t, x ; T, z) — (T — d g (x, z)) 2 — t 2 , y = (t, x) G M. x M n . (3.4) 

Note that as inj(M) > i 0 , (£, x) 1 —> ip(t,x',T, z) is smooth in M x (B g (z,i 0 ) \ {-})• 
Define 

S e , 7 = S(z,£,T,'y) (3.5) 

:= {(£, x) E [-T + £,T - £\ x R n ; if(t, x; T, z) > 7 2 , d g {x, z) < T}. 
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Also let z G M and 


E(z, £, T) = {(t, x) G M x M n ; |t| < T - £, \t\<T- d g (x, z)} (3.6) 

be the domain of influence of the cylinder 

W(z, £, T ) = (-T +£, T - £) x £ ff (z, £). (3.7) 

3.0.1 Some geometric estimates for domains of influences 



FIGURE 2. The hyperbolic surface between two domains of influence 

Lemma 3.2. LetT,£ : 7 be as in Assumption A5. Denote T^ n = ((T — I) 2 — y 2 )^ + 
Then 

E{z, e,T- 7 ) C S(z, £, T, 7 ) c E(z, £, T t „) U W(z, £, T) c £(*, 4 T). 

Proof. Assume that x G R” is such that s = d g (x,z) G [£, T — 7] and that 

\t\ < T — 7 — s. Then 

-0(L x; T, z) — (T — s ) 2 — f 2 

> (T — s ) 2 — (T — 7 — s ) 2 

> ((r — s) — (r — 7 — s)) ((t — s) + (r — 7 — s)) 

> 7 ( 2 (T — s) — 7 ) 

> 7 2 - 

We see that 

S(z,£,T, 7) D {(f,x) G [-T + T - £} x R n ; d g (x, z) > I, \T — 7 — d g (x, z)\ > |t|}. 
This proves 

E(^,T-7)CS(M,T,7). (3.8) 
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Assume next that s = d g (x, z ) G [£,Tg n \ and \t\ > Tn n — s. Then T ^ 7 — T < 0 implies 
that 

il)(t,x;T, z) = (T — s) 2 — t 2 

< ( T _£_( s _£))2_ ((rA7 _ £) _ (s _ £)) 2 

< ((T - £) 2 - (T £i7 - £) 2 ) - 2(T - f)(s - £) + 2(T £ , 7 - £)(s - £) 

< 7 2 - 

Thus, we see that the complement of S^ 7 satisfies 

{(t,x) G [-T + f,T-f]xr ; \T^-d g (x,z)\ < |t|} C S!„\JW{z,£,T). 

Hence we see that 

\ W(z, £, T) c E (z, £, T* i7 ) \ W(z, £, T). 

This and (13.8[) yield the claim. 

□ 


3.1 Applications: Stability on the domain of influence of a 
cylinder 

Here we consider the case when the solution of the wave equation ( 11 . 11 ) vanishes in 
the cylinder W{zq,£,T) and T may be so large that we have to consider also singular 
points for d g . We refer to Definition 13.II for the definition of sets used. 

Our aim is to prove the following: 

Theorem 3.3. Under the conditions of Assumption A5, let zq G K n , and define 
Q = (-T,T)xM, T = W(zo,T ,£), A = S(z 0 ,£,T, 1 )\T, 

D 0 = S{z 0 ,£,T, -J=) \ {(t,x) : t G M,d g (z 0 ,x) < ^}, Oi = O 0 \ T. 

Assume that u G i/ 1 (D) satisfies 

P(x,D)u(y) = f(y), foryefl 

and 


u\w(zo,e,T) ~ 0 . 


(3.9) 


Then for every 0 < 9 < 1 we have 

|MU 2 (A) < Ci63 




(In (e + Hlgl(ni) b e 
V m V e+ ||/lh 2(S2l , )) 


Here, Ci 63 depends only on a 0 , b 0 , b 3 , T, 7 , £, i 0 , and 6. 


Corollary 3.4. By Lemma ldTB we observe that, after a reparametrization of the time, 
E (zq,£,T) C S(z 0 ,£,T + 7,7). Consider the wave equation formulated in Theorem 
\d.M Hence for each 7 such that 0 < 7 < T — £, the optimal time of the control T — £ 
(with T — £ = raax Xty< =x( Zo / t T)\w(zo/,T)d g (y,x)) can be approximated from above by 
T — £ + 7 , using a result of stability of the unique continuation. 
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3.1.1 Local stability estimate 


Below, we say that the cut-off function corresponding to a center point y = (£, x) G 
Rxl” and a radius r is the product of a “time-variable cut-off function” and “space- 
variable cut-off function”, given by 




(3.10) 


<$(<)=/ l,i) (LA /# ) w=/ ,i) ( £ p £ 


), 


where G C*o(M n ) and G We assume that 0 < /^ < 1 and 

0 < f( tl ) < l. We assume that 


supp () C B r (0, a/2), / (a) |s K (o,i) = 1, 

supp (/ (sp) ) C B Rn (0, a/2), / (sp W(o,i) = 1. 

Then we have for y G M” +1 and r > 0 


supp {by.f) C C(y,2r), = 1- (3.11) 

Note that by (13. ip and (13.31) . C(y, 2 r) C C g (y, 2 y/btfr) and C g (y, -y/a^r) C C(y,r). 


For the proof of the global stability we must define the following points y and func¬ 
tions ij}~ f. 

Definition 3.5. (see Figure 2 below, where y = y r ) Let y = (x,t) G 

T, 7 )\{y; (GR; d g (z 0 , x) < £\. 

a) If t < d g (z 0 ,x ) < |i 0 , then define z = z 0 , T = T, ^f(y) = ij){y ; Zq,T), 

b) If |*o < d g (z 0 , x), then let (y) = ^(y; 'z, T ). Calling 7 ^ ^-a distance minimizing, 
unit speed geodesic from z 0 to x in M, we define z, T as follows: 

L:=d s (z„,x)-'‘j> 5 -^, T :— T — L, ?:=7 Wl f(2). 

Note that the choice of the point 'z is not unique as there may be several distance 
minimizing geodesics from z 0 to x. 

Observe that for £ < * 0 /4 and T > we have 

7 + «<7 +j<f = r- <4(zo, x) + | < T - ^ < T - l - (3.12) 

We then introduce the sets: 

Q 2 (z,T,£, 7 ) = S{z,T,£,^) D {(t,x) : t G M,£ < d g (z,x) < ^* 0 }, (3.13) 

O 

tt 3 (z,T,£, 7 ) = S(z,T,£, -^=)n{(t,x) : t G R, ^ < d 9 (?,x) < ^o}- 
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Lemma 3.6. Under the Assumptions A2, A5, let y and b e as Definition \3.5[ 
Then there exist r, R, cjj, cp, ca, and Ci 50 , C 131 , Ci 32; dependent only on the parameters 
a 0 , b 0 , b 3 , T, 7 , £, and i 0 such that the following property holds: 

If v G and y > 1 is such that 

supp 0) n C(y, 2 R) C {ye C(y , 2 R); if^f(y) < ^ f (y)} (3.14) 

and 


\\ v \\H 1 (C(y,2R)) ^ C U, \\P v \\L 2 (C(y,2R)) < c Pi II A(D 0 /fl) (bfrpPv) || L 2 ( R n+i) < Ca& ^ , 

then, 

2 

II-4(A)/ uf){by, r v) ||/fi(R«+i) < Ci 50 exp(-ci32/i“ ), for all u < y a /(3c 131 ). 

Proof. In the Appendix we have calculated uniform estimates for the function fi(y, Zq, T) 
defined in Definition 13.51 a): see ( 14.61) . (14.71) (with 7 j = 7 /\/2), (14. 8j) . (14. 121) . (14. 131) . 
Analogously, one can estimate the functions fi>(y,'z,T) defined in Definition 13.51 b): 
calling fl 2 = D 2 {z, T, £, 7 ) and D 3 = D 3 (^, T, £, 7 ), we have 


llc°(o 3 ) < h(T + 1), 

HV’-f||c°(n 3 ) - ^{T + l)((^/4) 1 + 1), 

11^1100(03) < W(T + l)((£/4)- 2 + m- 1 ), 

min \dij> $ f{y)\ > V^'fbfi 1/2 , min | p(y, dip~ f (y))\ > 2 7 2 , 

yeii3 ’ 3 

dRx(R n e)(^ 2 , dD 3 ) > —= mini — , ——-——, — 

1 ,en ’ Vbo 1 4 8 (T — £) 4 J 


(3.15) 


Next define i ? 0 = (2-sAo ) -1 min{i 0 /4, 3£/4, 7 2 /( 8 (T — £))} a uniform radius that let 
the ball T> 2 _r 0 (?/) stay inside the injectivity radius (in order to assure the regularity 
of if-zf) and inside the set D 3 (in order to assure that is non-characteristic in 
the ball), according to Lemma 14.31 Moreover, C(y,2R 0 ) C D 3 . We then consider 
the procedure of Appendix A to determine the radii r, R related to the function 
= , if~f(y)—'ijj~p(y). We set R\ = min{l, R 0 , (Alc^lc 0 ^)) -1 } hi the Table (14.31) . and 
we observe that, using the estimates (I3.15P for the derivatives of and Assumption 
A5, we can choose radii r, R , R 2 that are the same for each y , and consequently also 
the derived parameters. 

As seen in section^ all the parameters in Lemma 12761 r, R , cp, Cp, ca, < 750 , Ci 3 i, Ci 32 , 
depend on the uniform estimates for the quantities listed in Proposition 12.41 As we 
saw above, these estimates depend on the parameters ao, b 0 , b 3 , T, 7 , £, and Then, 
for each y , the claim follows from Lemma 12.61 for v in place of u , with the function 


3.1.2 Global stability estimate 

Rule of choosing the center points of small balls: 

We are going to apply the local stability estimate for the solution u of the wave 
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equation. Let r, R be the radii defined in Lemma 13.61 and consider the cylinders 
C(yj,r/ 2) having center points at y 3 chosen iteratively below, see (13.3[) . For each 
point y 3 we define a localizer bj(y) associated with y 3 = (t 3 , x 3 ) (see (13.111) ') by 

bj(t,x) = by jjr / 2 (t,x). 

We proceed in analogy with Assumption A4, with A in place of Q a , with A and fio 
defined in the statement of Theorem 13.31 Next set tp(y) = il>(y;T, z 0 ), as in (13.41) . 
The main difference is that here ip is not everywhere a C 2 ’ p function, as explained 
below. 

We define the set £ = {y 3 G A, j = 1,..., J 0 } and the functions Uj(y), iteratively as 
follows: 

1) For j = 1 we define u\{t,x) = u{t,x ) and consider A C O 0 

Vi e argmax {ij}(y; T, z 0 ); y = (t, x ) G A}. (3.16) 

2) For j > 2, after y±, y 2 ,..., y 3 -\ G £ have been chosen and the function Uj(t, x) has 
been constructed we proceed as follows: If supp (u 3 ) D A ^ 0, we choose y 3 to be a 
point that satisfies 

y 3 G argmax T, z 0 ); y = (t, x) G A\ C(y k , r/2)} (3.17) 


and define 


My) = (1 - bj-i(v))Uj-i(y). 


(3.18) 


We notice that by construction supp (u 3 ) DC(yj, 2 R) C {y\ ij}(y, T, z 0 ) < ijj(yj', T, ^o)}- 
When supp (u 3 ) D A = 0, we end the iteration and we set J 0 equal to j. 

Next we estimate the number of iteration steps Jo- By construction, the points y 3 in 
steps 1 and 2 satisfy d Rx (Rn, e )(yj, Vk) > r/2, j ± k. Moreover, 

— 2 

AcC( 2 //,po), where p 0 = —=(T + 1), y r = (0, z Q ), 

\/ai 

see (I3.3[) . Thus the maximal number Jo of steps is smaller or equal to the maximal 
number of points in a r net in the set A that is bounded by 


j < voIrxB " 1 (B(0, Pq + r)) (T + 2 ) n+1 
°~ vo1k xR »(C( 0,J)) _ 1 r n+1 


(3.19) 


where C\ is a uniform constant that can be estimated in an explicit way. 

Note that above we have always chosen y 3 = {t 3 ,x 3 ) as maximal points for the 
hyperbolic function i/;(y,T, zq) associated with the “original” center point zq and 
time T. The motivation for this choice is that the level sets of ijj(y, T , zq) are the best 
approximation of the domain of influence T.(zo,£,T) that we want to approach. 
When the distance of x 3 to the point zq is larger than the injectivity radius, the 
function y i->- ip{y, T , z 0 ) is only Lipschitz-smooth but it may happen that it is not C 2 - 
smooth. To apply Lemma 13.61 in this case, we choose a different hyperbolic function 
that changes at each step of the iteration and depends on the point y 3 . 
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FIGURE 3. Preparations to do unique continuation near yj = ( tj,Xj) G M x M n . 
When dj = d g (zo,Xj ) is larger than the injectivity radius, the boundary of the ball 
B g (zo,dj) may be non-smooth (the black external contour in the figure). For Xj G 
B g (zo,dj) we choose some distance minimizing geodesic 7z O) £,-([0, dj]) that connects zo 
to Xj. In the figure, this geodesic is the red curve from zq to x 3 . On this geodesic we 
choose a point Zj = 7 z 0 £j(dj — s o). The boundary of the ball B g (zj, so) (the red circle 
in the figure) is smooth and contains the point Xj. We do unique continuation near 
the point y 3 using the hyperbolic function associated with the center z 3 , that is 

smooth near y 3 . 

We then distinguish two cases as in Definition 13.51 

a) If I < d g (z 0 , Xj) < If-, then we consider if(y, Tj, Zj ) with Zj = Zo and Tj = T. 

b) Next, assume that dj := d g (z 0 ,Xj) > ^f.. Then, we define y = yj and as in 
Definition 13.51 

L i = L, Tj = T , Zj = z, Tj , zfi). 


Note that the choice of the point Zj is not unique as there may be several distance 
minimizing geodesics from zq to Xj. 


Lemma 3.7. For the points yj = ( tj,Xj ) G A, Zj G M n , the time Tj > 0, and 
the function Uj chosen above, the support condition \3.14\) is valid in the cylinder 
C(yj, 2 R) for the function ip ZjtTj (y), that is, 


supp 


{u j )GC(y jl 2R) C {ye M n+1 ; ^ Zj , Tj (y) < Zj )Tj (yj)}. 


(3.20) 


Moreover, we have fi>(yj', T, z 0 ) = fi>{yj] Tj , zf) and 


yj G dS{zj,Tj, i, 7 j) n C(yj, 2 R), where 7 := > 7 . 

Proof. If I < d g (zo,Xj ) < Ao, then the property is trivial because of fi(y, T)-, Zj) = 
ip(y; T, ^ 0 ) and the definition of yj G £ C A. 

Assume now that d g (z 0 ,Xj ) > Recall that by dehnition of R < R 0 in the proof 

of Lemma [3761 we have 2 R < b 0 ^ 2 < K Let us consider (t,x) G C(yj,2R). 
By the triangle inequality and the dehnition of Lj we have 


d g (x, Zj) + Lj > d g (x, z 0 ), d g (xj, Zj) + Lj = d g (xj, z 0 ). 
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This yields 


T - (d,(x, Zj) + Lj) < T - d a (x, z„). 


(3.21) 


Since 


d g (x, zj) + Lj < d g (x, Xj ) + d g (xj, Zj) + Lj < d g (x, Xj) + d g (xj, z 0 ) 

< 2Rb^ 2 + T — 7 < T, 

we have T — (d g (x, Zj) + Lj) > 0, and (13.2ip implies 

x] Tj, Zj) — (T — Lj -d g (x,Zj)) 2 -t 2 < (T — d g (x, z 0 )) 2 -t 2 = ip(t,x;T,zo).( 3.22) 
Hence, 

supp (■ Uj) fl C(yj, 2 R) C 

{(t, x) G C(yj, 2 R); ip(t, x; T, z 0 ) < 7 ?} C {(t, x) G C(yj, 2i2); ^(t, x; Tj, < 7 J} 

Moreover, when x G 7 ZOi ^.([L^, dj]), where 7 ^^. ([0, dj]) is a length minimizing min¬ 
imizing geodesic connecting zo to Xj, and i G K, we have d g (x,Zj ) + Lj = d g (x,zo) 
and 

il)(t,x;Tj, Zj) = (T — Lj — d g (x, Zj)) 2 — t 2 (3.23) 

= (T - d g (x,z 0 )) 2 -t 2 = ip(t,x’,T,z Q ). 

In particular, when (t,x) is equal to y 3 = ( tj,Xj ), we see that jj(t 3 , x 3 ; T 3 , Zj) = 
ip(tj,Xj-T,z 0 ). 

The above implies that 

C{y j ,2R)nS(zj,Tj,£,'y j ) C C{y j: 2R) fl S(z Q , T,i, 7 ,-). (3.24) 

Note that the boundary of S(zo,T,£,'yj) may be non-smooth in the ball C(jjj,2R), 
while the boundary of S(zj,Tj,£,'yj) is smooth. That is why we have introduced the 
new function ip Zh T 3 - We also recall that C(yj,2R) C S(z 0 ,T,£, -J=) and C{jjj,2R) C 
s(z i ,T j ,(,f 2 y 

By the construction of Uj and its support and the inclusion (13.24)1 we deduce that 

Uj = 0, for y G C(yj, 2R) fl S(zj, Tj, i, 7 j). (3.25) 

Moreover, since ip(tj, Xj\ Tj, Zj) = i/)(tj, Xj\ T, z 0 ) = 7 ?, we have yj G dS(zj, Tj, i, 7 j) fl 
C( yj ,2R). □ 

Proof of Theorem 13.31 We apply Theorem 11.21 in a special way. As mentioned 
before, here Lemma [3.61 replaces Lemma [2.61 
Set yj like in (13.171) . and u 3 like in (13.18^ . 

Step 1. Within the injectivity radius. Let d g (zo,Xj) < 7io/8. Define like in (13.13P 
and Lemma [3.61 

^ 0,1 — ^ 3 (^ 0 ; T, £,')), 

^i(y) = i>(y,T,z 0 ), 

I Xj 

Ai = n 2 {z 0 ,T,£,^). 71 = T g- > 7- 
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Here VWx,i = (T — £) 2 , V ; rmn : 1 = 7 i- Observe that C n 2 (z 0 , T, i, 7 ), the set used 
in Lemma 13.61 to compute the uniform radius R. By construction every y £ C(yj, 2R) 
is such that | < d g (z 0 ,Xj ) < hence x —> d g (z 0 ,x ) is C* 3 —smooth in C(y 3 ,2R). 
Moreover 1/7 is O 3 (r2o,i) and hence regular enough to apply the local stability result 
of Lemma 13.61 Here we are in the case where y = y 3 is like Definition 13.51 a). 
The condition (13.141) is fulfilled by u 3 due to the initial assumption that u — 0 in 
W(zq, T) and the construction of u 3 step by step. Call N\ the number of points y 3 
used to cover Ai. If T < 7i 0 /8, then the procedure stops here. If also T < 5i 0 /8, 
then it is enough to use a fraction of io above to define fl 2 and Q 3 . 

Step j > N\. Case T > 7%q/8. 

Here we change Do ,j, j(y) and A j at each step. We have 2 cases: 

a) If y 3 £ Q 0 ,i \ A 1; then we simply consider Q 0j = D 0)1 and 'ijjj = jpi (that is C 3 

regular since d g (x 3 ,z 0 ) < 7i 0 /8). Here A j = VL 2 {zq, T, i, 7 ,-), with 7 ^ = y 3 ; T. z 0 ), 
but Sj = {yj}, in the sense that we apply the local unique continuation step just 
once, in a cylinder C(y 3 ,2R) centred in y 3 £ {y;Vh(y) = (Uj)}- Observe that 

A j C O 2 (^o, T, l, 7 ). Again for y = y 3 holds Definition 13.51 ai and the condition (13.14)1 
is fulhlled by construction. By Remark 12.81 2.b) there is no need of defining i}>min,j 
here. 

b) If y 3 <£ D 0 .i- Then, d g (z 0 ,x 3 ) > 7i 0 /8. 

Here we are outside of the domain where "01 is certainly smooth, since the function 
x —> d g (z 0 ,x ) can fail to be C 2 —smooth in C{y 3 ,2R). So even if y 3 £ {y;Vh(y) = 
ipi(yj)}, to apply the local stability we choose another function ip 3 passing through y 3 
and having the good properties outlined in Lemma 13.71 Calling y = y 3 and defining 
Zj,Tj,^jj as in Definition 13.51 b). we can consider 

^ 0 ,j — ^3 { z j j Tj , £■, 7 ) j 

^j(y) = ^(y, T j,z 3 ), 

A j = 0 2 (z j ,T j ,^, 7 j ), 7 j = yj ip(yj;T, z 0 ) > 7 . 

Observe that A 3 C fl 2 (z 3 ,T 3 ,£, 7 ). Again we have S 3 = {yj}, in the sense that we 
apply unique continuation just in a cylinder C(y 3 ,2R) centred in y 3 £ {y; 1/7 (y) = 
i/)j{yj)}- The condition supp (u 3 ) nC(y 3 ,2R) C {y;$/(y) < i’ji.Vj)} is fulhlled due to 
Lemma [3.71 By Remark 12.81 2.b) there is no need of defining VWiy here. 

Notice that, due to our uniform estimates, the radii R and r remain unchanged for 
every y 3 and the other constants of the Table (14.3(1 are chosen uniformly. This implies 
that C 156 ,* = Ci 56 ,i- We also recall that C{y 3 ,2R) C S(zo,T, i, -^=) for every j , by the 
construction of the points y 3 and the choice of R. Here Z(y) £ C'^°(M n+1 ) satishes 
/ = 1 on U J k ° =l C{jjk,2R), 0 < l < 1 and supp(Z) C S(zo,T,£,-J=), see Remark 12.81 4. 
The coefficient C 163 is computed like Ci 6 i- □ 

Remark 3.8. We remark that an alternative proof of Tli. 13.31 is possible by applying 
Th. 11.21 in the following way. Define a net of center points {t k ,z k ) for the translated 
hyperbolic functions: 

V>(y; T k , z k , tk) = (Tfc - dg{x, z k )) 2 - (t- tkf. 
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Choose T = W(z 0 ,,T,£), O 0 ,fc C {y,y E [~T k + t k ,T k + t k ] x M n ; i()(y,T k ,z k ,t k ) > 
7 fc/ 2 , T k > d g (x,z k )} and A k C {y\ y G [~T k + t k ,T k + t k ] x M n ; if>(y,T k , z k ,t k ) > 
T k > d g (x,z k )}. The construction is similar to the one in Figure 1 of section 
[2j In this case one does not need to introduce the points y of Definition 13.51 The 
parameters (t k , z k , T k , 7 *,) should be chosen such that Do,fc is contained in the domain 
0 < d g (z k ,x ) < (to guarantee the regularity of if(y;T k , z k ,t k )). Moreover A k C 
T,(z 0 ,£,T) and their union should cover a subset of the domain of influence T,(zq,£,T). 


3.1.3 The case of solutions with small values in a cylinder 

Our purpose is to reformulate Theorem 13.31 for a wave equation with vanishing source 
term and a solution u that is no longer vanishing but it is small inside a cylindrical 
set. 


Corollary 3.9. Under the conditions of Assumption A5, let Zo G M n . Also, let 
n = (—T, T) x B g (z 0 ,T + £), n, = S(z 0 ,£,T,^)\{(t,x) : t G M,d 3 (^ 0 ,x) < if A}, 
O 2 = S(zo,£,T,^f). Assume that u G H l (Cl) satisfies 


P(x, D)w(t, x) = 0, for (t,x) G Cl 

and define 

W\ — (—T + £,T — £) x B g (z 0 , £ + j). 
Then for every 0 < 9 < 1 we have 

||w||fl-i(ni\w 0 ) 


\ w \\l 2 (C1 2 \Wi) < c 166 


(ln(e+ )) e 

V V C'\\w\w l lOpM/p'' 


(3.26) 


(3.27) 


Here, Ci 66 depends only on ao, b 0l b 3l T, 7 , £, i 0 , and 6. 

Proof Let B 0 = B g (z 0 ,£), B\ = B g (z 0 ,£ + 7 ), and W 0 = W{z 0 ,£,T) = (-T + £,T- 
£) x B 0 C W 1 . We use a cut-off function rj{x) G Cq(B 1 ), 0 < 77 < 1, that is equal 
one in B 0 and satisfies ||t7||c 2 (R") < c o7~ 2 , where c 0 is a uniform constant. Then 
w(x,t ) = (1 — rj(x))w(x,t) vanishes in W 0 and we have 

P(y,D)w(t,x) = F(t,x ), in 0, 

F{t,x) = —g jk (x)(DjT]D k r])w — g jk (x)(DjgD k w) — h j (x)(Djg)w G L 2 (0), 
and since 77 is supported in B\, 

||-F||l 2 (Oi) < ||-F1|z, 2 (Wi) < ciUwlwill^qwi), (3.28) 

where C\ is a uniform constant. Also, since w — w in 0 \ W\, we have 

||w ~ Hl-ffhW) — c 2||' u; |w r i||rf 1 (Wi)) 

IMI-ffhfii) < c 2 ||H|.ffi(n 1 ) 
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(3.29) 

(3.30) 








where C 2 is a uniform constant. Summarizing, above we have seen that 


P(y,D)w(t,x ) = F(t,x), in (-T + £,T — £) x B g (z 0 ,T + £), 
w\ Wo = 0 . 

Moreover F in (— T + £, T — £) x B g (zo,T + £) vanishes outside of W\ and F is small 
if ||^|TViH-ffhWi) is small. Then, applying Theorem 13.31 to w we get 


MU 2 (RAWo) < c 163 


\w 


ll/ZpOi 


\Wo) 




H F Hi 2 (n 1 \iv 0 ) - 

As \\w\\ L 2 {Q2 \ Wl) < \\w\\L\n 2 \w 0 ), and by (ESQ}, 

IMItfRfiAWo) < c 2 ||u'|| H i(n 1 \w 0 ) ! 

and since the function t is increasing for t > 0 , we get 

C2||w||//i(Q 1 \vV 0 ) 


MU 2 (H 2 \Wi) < c 163 


(In (e + C 2 iM llH M (n i xw °) )) 9 ' 

V V Cllliulwi II zrl/ixr.w / / 


c ill ul lu , 1 lliri(v^ 1 )) 

This proves the claim with Ci 66 = max(c 2 /ci, C 2 C 163 ). □ 
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4 Appendix 

4.1 A: Geometric constants 

We write the table of the constants used in the article. This is a special version of 
Table (5.1) in |j], since now all the coefficients are calculated independently of the 
center point r/k and of the local information. 

In order to get the uniform coefficients we use the same notations as in Section 3.1 
of ( 6 ]: 

a) By Assumption Al, we consider the case of the wave operator (II. ip with princi¬ 
pal symbol p(y, f) = -£o + Ej fc =i 9 jk (x)tj£ k , with 0 < a x 5 jk < g jk (x) < 6 , 6 jk , 
ai, b\ > 0. 

Call e = (£ 0 ,0 e R x R n , where |£| 2 = £” = i £?. 
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b) We consider the function 0 G C' 2,p (M n+1 ), for some p G (0,1), such that 
pil/i'ip'iy)) 0 0 and 0'( y ) ^ 0 in a domain f2 0 C fl. Let y 0 G h2 0 be a general 
point lying on the level set S = {y\ i^(y) = 0}. Call pi = min yg ^ o p{y, ip') > 
0, Ci = min^ o \ip'(y)\>0. 

Moreover we use Einstein’s convention for the repeated indexes. 

We recall the three Steps - procedure to calculate the geometric parameters in [6). 

Step 1 . Given a function 0 G C' 2,p (R n+1 ) fulfilling the assumptions above in a domain 
we find positive constants M 2 , Mi, Mp such that the following inequality 
holds true 


+ Mi ( \ p ( y ’ + \{p'i{y^ + ir^\y),^\y))\ 2 ) 


{p{y,€ + i ^ , (y)),p(y,t + iT'ilj , (y))} 

2 IT 


> Mp(r 2 + |£| 2 ) 


for every £ G R x R”, ( / 0, t G 1. The previous inequality proves that the 
hypersurface S = {y; ip(y) = 0} is conormally strongly pseudoconvex w.r.t. P 
in 


Step 2 . For 0 = e A 0 with y 0 on the level set 0(?/) = 1, we find A > 0 such that the 
following inequality holds true 


M 2 £q + 


Mi 


I p(y,t + iT(jf(y))\‘ 


min{l, A 2 4> 2 (y)} 


t 2 + |£| 2 

1 {p(y,Z + iT(f/(y)),p(y,£ + iT<j/(y))} 


M(y) 


2 ir 


> M P min{l, A 0 (y)}(r + |£| 


for every £ G R x R n , £ / 0, t G 1. The previous inequality proves that the 
function 0 is conormally strongly pseudoconvex w.r.t. P in f2 0 . 


Step 3 . We consider a perturbation of 0 by the shifted 2nd order polynomial centred 
in the point y 0 , 

f(y) = (-^o) (y - yoT/v\ - a\y - y 0 1 2 . (4.1) 

M<2 

In a ball B(y 0 ,Ri) C O 0 where /' 0 0 we define 

0o = min 0(?/), 0 m = max <f>(y). 
y&B(y 0 ,Ri) V&B(y 0 ,Ri) 

We hnd a and R 2 > 0 small enough such that in the ball B(y 0 , R 2 ) the following 
inequalities hold true: f(y) < 0(y) in B(y 0 , R 2 )\{y 0 }, and 

m <-2 , o \p(y,Z + irf'(y ))| 2 , {p(y^ + iTf'(y)),p{y^ + iTf (y))} 

+ 2Ml —— + -- 

> - 2 (r 2 + Id 2 ) . 

The previous inequality proves that the function / is conormally strongly pseu¬ 
doconvex w.r.t. P in B(y 0 ,R 2 ). 
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How to obtain the uniform Table (14.3f) . starting from the table in [6] 

We need to recalculate i? 2 invariantly. Hence we take: 

a = 2 c t R p 2 , c t = nlA^Ua^.Oo (4.2) 

|A^| max,n 0 = 4 >m max(A|'0 // |co,p ( n o) , A 2 |^| c o,i(n 0 )|V / '|c°(n 0 b ^Mc^nojIV’Icopio)) 

and after a first estimate in B Rl (yo) we estimate Mr calculating the norms in Ho 
instead of in B Rl (y 0 )\ 

Mr > 1 — 

+ |AV J |mo i c,no-^2^^(l + ^VmIV’IcoIHo) + W C°(O 0 ) + ^ 1^IC°(t2 0 )))] ■ 

Recall also that 


\f'\c°(B R2 ( yo )) < A0 mW\c°{q. 0 ) + ^n\\ijj\ max> Q. 0 R^ rP , 

\f”\c°(B R2 (y 0 )) < A^MdV'^lc 0 + -^|'0 / |c , °) + 4n|A'0|max,n o -^2) 

l^(v)l = l<M« + AV/ X ^)l > e _1 A 2 Cf/2. 

In the table we anticipate the dehnition of eo in order to embody in R 2 the condition 
R 2 < (16 + i^) 1//2 ; that is essential to define the inequality in Th. 12.31 The es¬ 

timate for the minimum r was done in [6] and we refine it calculating the norms in Ho- 

Consider now the points yk G £ defined in Assumption A4. For each k , via the 
translation ijj(y ) — ijj(yk), we have ^(Uk) = 0 and we can replace yo with yk in the 
definitions above. After the previous translation, i.e. cj)(y) = exp(A {^(y) — ijj(yk))), 
we still have e" 1 = </>o < (j>(y) < 4 >m = e in each ball B Rl (yk) and we can consider 
<f>(yk) = 1- All these translated curves share the same parameters as in Table (14.31) . 

Remark 4.1. In [6] we assumed that distRn+i{5H 0 , <9fi} > 0. 

Here we have the condition dist®n+i{(9Ho, fi a } > 0, and then we assume that C flj. 
In Section [3] we apply Table 14.31 to a C 3 — function A- Instead of calculating the 
C 2,p — norm, it is more practical to use the C 3 — norm, but this requires the following 
modifications. Set p — 1 and c R in place of n\\i/j\ max tf 0 , where 

W\c°{Br 2 ) < Cti ■= Mm{W\c°(P.o) + 'M'0'|co(n o) ), 

W\C°{Br 2 ) < C T 2 := A0 M (3A|^ / |c°(Oo)lV ;// |c o (no) + -^ 2 |^ / |cO(Oo) + W\c 0 (p. 0 ))i 

R 3 1 1 

cr = 2 ctR 2 , ct = cti + c R 2 , 5 = c R q 2 —y, q = -(16 + —) _1 ^ 2 , 

8 4 16 

f — 4> < —c T q 2 R\, \(ft' — f\c°(B R2 ) < 5c t R%, \(p" — f" | c° (Br 2 ) < 5c t R2- 
The bound for r remains unchanged since IA'Ico^br )q 2 R%/8 < <5 (see |6j). 
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Table for the constants calculated as in ® (4.3) 


Name 


Limit Value 

c 3 

> 

20(1 + rrV% 1(f2o) ) 0' ci(n 0 )( 1 + W c°(n 0 )) 

M P 

< 

1 

M, 

> 

(M p + C 3 ) max{^, 2(p ' i)2 } 

m 2 

> 


A 

> 

max{Mi, e, ^ |^° (n o) j. 

00 

> 

e- 1 

4>m 

< 

e 

Ri 

< 

mm{l, dist R »+i{dfl 0 , ft a }, A|^lco ( n n) ^ 

Cioo 

> 

10(1 + 

eo 

< 

I 


R-2 

< 

min Sr, ( _ Cl _ ) ( x2 * mC ? )p 

l X ’ \ 20 M (|l/'"| C O( r!o) +A|V’'lpO ( Q o )+ 4n l Al /’|ma a; .n o /A)/ ’ V4n|AV>|max,n 0 ' ’ 

( 1 A 2T2 p 1 eo (1 « 1 l \ 1 / 2 

V4doo|A^IL^ M X (1+A^IV|* 0(no) ) ' ’ 8 ^ l iU ^ 16 ) ’ 

( A?io A h 

4cioo|A^| ma:Ei n 0 (l+A 2 </>| f |i/i , |^, 0 ^ nn ^+A 2 (/>| f (|i/) , | c , 0 (n 0 )l'i/ , "l c . 0 (n 0 )+A|'i/i , |^ ;0 ^ rin ^) 

a 

> 

277-1| rnax,£l o -^2 

T 0 

> 

ma,x{l, 64(4Mi + 4A ^J ((A0 m(|0"|c° + A|0'|^ o ) + 
4n A0 max&o -^2) 2 (1 + ^ 1^ 1 C°(O 0 ) Y + ™ 1' h 1 C0(Q 0 ) ( 2 + 
2(A0 M |0 / |co(n o ) + ^n\\ip\ max ^ 0 R\ +p f + 2|g|^ 0(no) )} 

R 

< 

ki6+A)' 1/2 ^~ 

6 

< 

n^(l6+^) 1 A0| ma;E) no-R2 +P 

r 

< 

nA a Cf J(l6+^) _1 J^+'’ 

2e (| </>'| c o (n n ) +5n| 0" | c o,p (f2n )) 

Cl,T 

> 

\/ 4 ( r 0 + 4(A</>o) ) 

C2,T 

> 


Cl33 

> 

2(1 + n \g j |c°(n 0 ))( ro ( 4 i?) 2 1 m (-*- + V' c»(fi 0 ) + 

5n|A0| ma ^ o R 2 +P + ^ m )) 


4.1.1 Uniform estimates for the hyperbolic function ip 

Uniform regularity estimates for the distance function d g . It is a well known 
fact, see [9], that if a metric is C' m -srnooth, then the Riemannian normal coordinates 
are C ,m_1 -smooth, and the metric tensor in these coordinates is C ,m ~ 2 -srnooth. In 
particular, the distance function x i —> d g (x,z) is C' m ^ 1 -smooth. In the following 
we consider how to obtain uniform bounds for the distance function under suitable 
assumptions. 
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Let m > 2 be an integer, and a, r 0 > 0 , Q 0 > 0 be fixed parameters. 

Assume that on M are local coordinates (Uk, '&&), where Uk C M are open and 
T k : Uk —> M n such that: 

(i) For any x G M there is k such that the metric balls B g (x,r 0 ) C Uk, Let Wk = 
*k(Uk) C M n . 

(ii) The metric tensor satishes in local coordinates 

J/< (**).<?< 4/, (4.4) 

(iii) We have ||(’L fe )*e/|| Cm(Wfc) < Qo, 

(iv) The transition functions satisfy: ||4/ fc o < Q 0 , 

(v) The injectivity radius satishes: inj (M,g) > 2r 0 = Iq, where 0 < < oJkTi anc ^ 

A m is an upper bound for the sectional curvature of M. 

Under the previous assumptions and using the notation (v) = (1 + Ini 2 ) 1 / 2 , one 
can obtain the estimate for the derivatives of d g when 0 < d g (x o,x) < i o: 

\D°d g (x 0 ,x) \ < e A»Wo> 4 ftN,4> 2 ( | s ( X o, x )( 1 -l“l), | a | > 0. (4.5) 

Here c^ n > 1 are coefficients which depend only on m, n; their value may be explicitly 
found from combinatorics. 

Consequently, we consider Assumption A5. Then for all z G M n and A\ = {x G 
IR n ; \t < d g (x,z) < |*o} we have (see (14. 5 j) and [9] for details) 

\\d g {-,z)\\ c i(A 1 ) <b 2 , \\d g { ■ ,z)||c 3 (Ai) < b 2 (£/A)~ 2 , (4.6) 

where 62 depends on ao,b 0 ,b 3 , and io in an explicit way. 

Let z G M, T > 0 and recall the ’hyperbolic function’ introduced in Definition 13.11 

■0 g ,T(t , x) \= ip(t, x; T, z) — (T — d g (x, z)) 2 - t 2 . 

In the following we consider properties of this function in order to construct the 
related Table (14.31) . 

We recall that the principal symbol of the wave operator P, at y = (t, x), £ = 
(£ 0 , £ 1 , • • •, 60 e T*(R x W 1 ) by p(y, £) = -f 2 + ^= 1 9 jk (x)Uk- 

Lemma 4.2. Let z G K n ; A\ = {x G j^£ < d g (x, z) < |i 0 } and A = [— T, T]xA 3 . 
Let y = (t,x) G A be such that |f| < T — d(x,z). Also, assume that ip Zt T{y) > 7 /, 
with 0 < 7 1 < T. Then d'lpz^iv) Is well defined and satisfies 

p(y,dij> z ,T(y)) > 47 ?, (4-7) 

min \dfi Z:T (y)\ > 2'j I b^ 1/2 . 

y&A 

Moreover, we have 

II^2,tIIc 0 (.4) < h(T + 1), 

Wz,t\\c°(A) < bfiT + l)((^/4)" 1 + 1), 

Kt\\c°(A) < h(T + l)((£/4)- 2 + {1/ 4)- 1 ), (4.8) 

where 64 = 64 ( 00 , 60 , 63 ) is a uniform constant. 
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Proof. Consider the co-normal of the level set of if Zj T, 

v = di/j ZtT (y) = (-2 1 , -2 (T - d g (x, z))d x d g (x, z)) e T*(R x IT). 
Using the two following facts : 


y g 3k (x)d Xj d g (x, z)d Xk d g (x, z) 
j,k =1 


livd fl (-^)n; = i, 


- &+ y ^(x)^ < + y =■ yen 

j,k =i j,fc=i 


2 

dt 2 +g) 


we obtain 


P(?/, dif Z}T (y)) = 4(-t 2 + (T - d s (x, z)) 2 ) > 47 2 > 0, (4.9) 

|# z ,T(?/)|dt 2 +g > = 27 / implying |d^* )T (j/)| > 27 /&q 1/2 . 

Let us now analyze derivatives of if Zj T- We have 

VV7,T|(t,x) = (-2f, -2(T- d 9 (a:,2:))V x d 9 (a;,2:)), 

V^tIm = (-2, —2(T - d 9 (x,^))V 2 d s (x,^) + 2 \7 x d g (x,z) ® V^d^z, z)), 


V 3 -0z,T|(t,x) = (0,U), 

U = -2(T - d s (x, z))V 3 d 9 (a:, z) + 4V ;c d 3 (a:, z) ® V^d g (x, z) 

+2V x d g (x, z) ® V x d g (x, z ). 

Thus, by calling if' = df) z ^, if" = V 2 , 02 ,t, if'" = V 3 ifz,T, we get 

IIV ; 2 ,tIIc , 0 (^) - 4(7" + l)||d s ( •, ^)||c 1 (Ai), 

11 Virile 0 (. 4 ) < 4(T + 1)(11 dg (•, z) ||c 2 (Ai) + | \d g { •, 2)||ci(Ai))j 

Wz,t\\c«{A) < 6(1 + T)(||d 9 (-,^)|| C 3 (Al) + ||d s (-,z)|| C i ( A l) ||d s (-,z)||ca ( A l) )(4.10) 

and then one can use (14.61) at the right hand side, where 64 is a uniform constant. □ 

Next we estimate the distance between two level sets of if z ,T, or if~. f, outside of 
a cylinder of radius £ (see Definition 13.ip and its consequences. 

Lemma 4.3. Under the Assumption A5, we have, 

a) For i = 1,2 define Li — {y; y e dS(z, T, i, 7 *), d g (x, z) > ij, with 71 = 7 /-\/2 and 
72 = 7 . Hence, 

7 2 

dist dt 2 +g (L l ,L 2 ) > _ := Ci 80 . (4.11) 

Consequently, defining zq, A and D 0 as m Theorem, 1,9. .91 we get 

1 3£ 

dzsf R »+i(A, <9D 0 ) > - 7 = min{ci 80 , -r}- (4.12) 

V^o 4 
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b) Let y = (x, t) be a point defined in Definition ^. 51 and let fi(y, T, z ) be the associ¬ 
ated hyperbolic function. For i = 3,4 define L * = {y,y £ dS(z, T, £, jf), d g (x, z) > £}. 
Then, for 73 = 7 / 1/2 and 74 = 7 , 


dist^+ g (L3, Lfi) > ci8o- 


(4.13) 


Proof, a) Let yo = (to,xo) G L 2 and y\ = (ti,xi) G L\. Our aim is to get a positive 
lower bound for distd t 2 +s (?/ 1 , y 0 ). The two points lay on two level sets of fi>(y,T,z), 
and we consider their difference: 

'fiiVo; T, z) - T, z ) — (T — d g (z , x 0 )) 2 d g {z, x /) 2 + t\ = y 

{dg(z,x 1 ) - dg(z,x 0 ))(2T - dg(z,x 0 ) -dg(z,x 1 )) +t\ -t 2 0 = —. (4.14) 
By definition, we know that 


(2T - d ff (z, x 0 ) - ( 2 , xi)) > 0, £ < d g (z, x 0 ) < T - 7 , 


£ < d g (z, Xi) <T — 


7 

v/2 


< T. 


Assume w.r.o.g. that t 0l ti > 0. 
Case 1. f 2 — to = ^ [0,1), 


d fl (xi,x 0 ) > (d s (z,xi) - d fl (z,x 0 )) = 


(1 — m) 7 2 


(1 — m) 7 2 

2(2T - d fl (z, x 0 ) - d ff (z, xO) 2{2T-2£fi 


> 


Case 2 . tf — tl = y, that implies d g (z,x 0 ) = d g (z,x 1 ) and d s (xi,x 0 ) > 0. 
Case 3. t\ — t^ = q^-, q > 1 


(dg{z,x 1 ) - d g (z, x 0 ))(2T - dg(z,x 0 ) - d fl (z,xi)) = -(g - l)y 
here one can reverse the signs and prove as in case 1 : 

(? - 1 ) 7 2 


dg(x i,x 0 ) > (d g (z, x 0 ) - d 9 (z,xi)) = 


(g ~ !)7 2 

2(2T - d fl (z, x 0 ) - d g (z, xO) 2(2T - 2£)' 


> 


Case 4. t{ — t^ = —p\ < 0, for p > 0 


(d ff (z,xi) - d g (z, x 0 ))(2T - d g (z, x 0 ) - d fl (z,xi)) = 


(1 +p)7 2 


Case 1. is then the worse case. Hence for t\ = \Jm'Y^ + tf < \J-y + (T — £) 2 and 
since 7 < (T — £), we have 


(+ 2 _ 

dist fir ( 2 /i, 2 / 0 ) =max{|t 1 « f 0 | = 77 - 7 -^y, d g (xi, x 2 )} > 


> max{m 


7 


+ (T — £) 2 + T — £) 


, (1 -m) 


(ti + to) 

7 2 


}> 


7 


2(2T — 2 £) J - 8 (T — 
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b) We then consider yo,yi belonging to two level sets of the function if(y,z,T) to 
calculate the left hand site of (I4.13p . We repeat the same computation as above with 
the new values. We recall (13.121b By triangular inequality, since T = T — d g (z 0 , z), 

(2 T - d g (z , x 0 ) - d g (z, xi)) <2 T - d g (z 0 , x 0 ) - d g (z 0 , xi) < 2 (T - £). 

And again t x = + tg < + (T - f?) 2 < 2 (T - £) < 2 (T - £). 

Hence the new estimate is bounded from below by Ci 80 . □ 
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